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Preface 



This book is devoted to an investigation of the vacuum of quantum elec- 
trodynamics (QED), relying on the perturbative effective-action approach. If 
the vacuum is probed with external perturbations, the response of the system 
can be analyzed after averaging over the high-energy degrees of freedom. This 
results in an effective description of the properties of the vacuum, which are 
comparable to the properties of a classical medium. 

We concentrate primarily on the physics of slowly varying fields or soft 
photons by integrating out the high-energy degrees of freedom, i.e. the elec- 
trons, employing Schwinger’s proper-time method. We derive a new represen- 
tation of the one-loop photon polarization tensor, coupling to all orders to 
an arbitrary constant electromagnetic field, fully maintaining the dependence 
on the complete set of invariants. 

On the basis of effective Lagrangians, we derive the light cone condition 
for low-frequency photons propagating in strong fields. Our formalism can be 
extended to various external perturbations, such as temperature and Casimir 
situations. We give a proof of the “unified formula” for low-energy phenom- 
ena that describes the refractive indices of various perturbed quantum vacua. 
In the high-energy domain, we observe similarities between a vacuum with a 
superstrong magnetic field and a magnetized plasma. The question of mea- 
surability of the various effects is addressed; a violation of causality is not 
found. 

Furthermore, the QED one-loop effective Lagrangian in the presence of an 
arbitrary constant electromagnetic background field at finite temperature is 
derived. We show that the requirement for gauge invariance induces an addi- 
tional invariant on the quantum level, which is physically related to a chemical 
potential and to Debye screening. Further applications are presented. 

In our treatment of 2+1-dimensional QED, we discuss several unfamiliar 
features of the effective-action approach, such as the formation of a magnet- 
ically induced chiral condensate, the perturbative generation of a parity-odd 
Chern-Simons term and a derivative expansion of the QED 2+1 effective ac- 
tion for inhomogeneous fields. 

During the course of writing the manuscript, we have benefited from dis- 
cussions and correspondence with many colleagues. In particular, we wish 
to thank Professor Martin Reuter for numerous discussions and various use- 
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ful suggestions. We are indebted to Dr. Klaus Scharnhorst for his helpful 
comments, criticism and encouragement. 

We are also grateful to Mrs. Virginia Dittrich, who helped us to put this 
work into correct English. 



Tubingen, Walter Dittrich 

April 2000 Holger Gies 
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1. Introduction 



During the development of quantum field theory, it has become a popular 
point of view to consider the quantum vacuum as a medium. The manifold 
phenomena which arise from the presence (or the provoked absence) of virtual 
particles in the vacuum, such as the Lamb shift or the Casimir effect, tempt 
one to assign real properties to the vacuum. On the other hand, a physicist 
can hardly accept an establishment of this kind of “modern ether” without 
reservations; it is, at least, unsatisfactory to formulate a fundamental theory 
with the aid of some ingredients which elude direct measurement. 

In the present volume, we do not insist on either point of view as the 
one and only foundation of physics, but intend to demonstrate that it can be 
pragmatically appropriate and useful to understand certain quantum effects 
as peculiarities of the vacuum - both intuitively and formally. 

The formal concept of relating the full quantum theory to properties of the 
vacuum is given by the effective action. By integrating out the high-energy 
degrees of freedom of the exact theory, one arrives at an effective description 
of the low-energy degrees of freedom which are relevant to the physics of the 
vacuum. 

This procedure is particularly successful for theories such as quantum 
electrodynamics (QED), where it is easy to identify and separate the differ- 
ent degrees of freedom. In QED, it is the electron, which is separated from 
the low-energy photons by its mass, that represents the fundamental scale in 
the theory. In strongly coupled (gauge) systems such as quantum chromody- 
namics (QCD), one has to rely on some more or less justified approximations 
in order to establish the effective-action approach. 

The effective-action approach has also proved useful when additional en- 
ergy scales are introduced which can be varied over a wide range compared 
with the high-energy scale. These additional scales may represent pertur- 
bations of the vacuum, probing its response under external influences. The 
(pragmatic) idea then is to compare the response of the vacuum with the re- 
sponse of a medium under the influence of such a perturbation. The properties 
of a medium are finally assigned to the vacuum itself in order to complete 
the comparison. 

In the present work, we apply this philosophy to QED, which is pertur- 
batively accessible. We primarily employ external electromagnetic fields as 
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1 . Introduction 



the perturbation of the vacuum, i.e. as the additional energy scale. Further 
examples of perturbation are given by coupling the system to an external 
heat bath or limiting its spatial extension by means of perfectly conducting 
(Casimir) plates. Also imaginable are gravitational perturbations and per- 
turbations that arise from space-time embeddings with nontrivial topology. 

Technically speaking, we integrate out the fluctuating particles which are 
coupled to the external perturbation, to all orders if possible. The result is 
then compared with the unperturbed vacuum, and the effects that arise from 
the discrepancies are investigated. Since the remaining low-energy degrees of 
freedom of QED are electromagnetic fields, we study the properties of the 
modified vacua by, in particular, searching for a modification of the propa- 
gation of low-frequency light in these modified vacua. It is proposed that the 
deformation of the light cone can be used as a detector of the various char- 
acteristics of the vacuum. Our intention is that this theoretical work can be 
related to experiments, currently in preparation, that are aimed at proving di- 
rectly the quantum-induced nonlinear extensions of classical electrodynamics 
for the first time [167]. 

These investigations are presented from two different starting points. In 
the first part, we start with the complete quantum theory and study the 
polarization tensor and induced current in an electromagnetic field in detail, 
whereas in the second part, we take the effective action as the primary object. 
Once the effective action is deduced from the quantum theory, we can employ 
it as the definition of a new classical theory which effectively incorporates the 
quantum physics. 

In the perturbative approach, the proper-time method as developed by 
Schwinger is extremely useful, since it not only represents a powerful reg- 
ularization technique, but also supports a detailed investigation of gauge 
invariance. Moreover, it provides for a clear formulation and solution of the 
technical difficulties. 

Another important tool, which we use throughout these investigations of 
QED, derives from the fact that the effective Lagrangian is a scalar (in all 
respects) and can therefore depend only on the invariants of all allowed sym- 
metry transformations of the theory (gauge and Lorentz symmetry). Iden- 
tifying the complete set of invariants at the beginning of the investigations 
allows us to formulate the calculations in full generality, and simplifies the 
calculations in most cases. The simplifications arise from the fact that the 
computations of the dynamics of a system are separate from the algebraic 
computations referring to the symmetries of the theory. 



2. Nonlinear Electrodynamics: 
Quantum Theory 



2.1 Introduction to Proper-Time Methods 

We give a brief survey of the proper-time technique as formulated by Schwinger, 
in order to deal with quantized fermionic fields interacting with an external elec- 
tromagnetic field. The relation to one-particle quantum mechanics is emphasized, 
and questions of gauge invariance are studied. We describe three different meth- 
ods for computing the proper-time transition amplitude, which represents the cen- 
tral object of this technique. The derivation of the fermionic Green’s function and 
the Heisenberg-Euler effective Lagrangian of QED in Schwinger’s representation is 
sketched. 

The main purpose of this section is to establish conventions and notation, as 
well as to lay the foundations for subsequent investigations. 

One of the most important tools for the analysis of quantum electrodynam- 
ics is the proper-time method introduced by Fock [79] and Schwinger [148]. 
It has been proved extremely useful in studying propagators and effective 
Lagrangians, especially for the case of additional external fields. 

Schwinger’s main intention was to circumvent problems associated with 
gauge invariance. The proper-time method explicitly deals with objects that 
transform gauge-covariantly, thereby ensuring the invariance of the procedure 
and its results. In the literature, the proper-time method is in most cases 
simply used as a regularization procedure, because it is able to isolate the 
divergences of a calculation with respect to the proper time. Lorentz and 
gauge invariance of this isolation process is maintained, since the proper- 
time parameter is not related to a particular reference frame or choice of 
gauge. 

But more than being a regularization procedure, the proper-time tech- 
nique relates the field-theoretic problem of how a particle field interacts with 
an external (“c-number”) electromagnetic field to the description of a par- 
ticle’s motion with respect to an additional evolution parameter, the proper 
time (the “fifth parameter”). In this sense the technique falls back on the 
principles of classical and one-particle quantum mechanics. 

In the following, we shall give a brief introduction to the proper-time 
method, aiming especially at the construction of the Green’s function of a 
Dirac field and at the calculation of the effective Lagrangian. The electromag- 
netic field is considered to be an external one, i.e. it is treated as a c-number; 
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hence, radiative corrections are neglected and we are, by definition, dealing 
with the one- loop approximation. 

We begin by stating that the vacuum-to- vacuum persistence amplitude in 
the presence of an external electromagnetic field A M is related to the effective 
action (effective Lagrangian) by 

(0+ | 0_ ) A = e iW,(1) M = jfd'xcWlx)' (2.1) 



The superscript (1) indicates the one-loop character of the quantity to which 
it is attached. W ^ is defined so that it generates the vacuum expectation 
value of the current j M upon differentiating with respect to the external field 

A,: 



SW^[A] 

SA^x) 



= (0|j M (*)|0 Y 



( 2 . 2 ) 



According to Schwinger, the correct definition of the current operator is based 
on the explicitly symmetric treatment of the oppositely charged Dirac fields: 1 






(2.3) 



When the desired vacuum expectation value of the current operator is calcu- 
lated, the charge symmetrization translates into a time symmetrization, and 
we finally recover the definition of the propagator of a Dirac particle: 

(0|j M (a;) I0)" 4 = — e lim 7 M ( 0 | T |0) A 

x' — >X 

= ie trpy^ G(x, x\A)\ , (2.4) 



where the limit has to be taken symmetrically with respect to the time coor- 
dinate, and T denotes the time-ordering operator. The propagator G satisfies 
the Green’s function equation of a Dirac particle: 



Ih^n^) + m\ G(x, x'\ A) = S(x - x'), (2.5) 

with 77 n = —\d^ — eA^. Equation (2.5) can be interpreted as the rep- 
resentation of an operator equation in configuration space (G(x,x r |A) = 
(x\G[A]\x')): 

(yTI + m) G[A\ = 1. (2.6) 



We introduce the proper-time representation of the formal solution of the 
operator equation (2.6), 



OO 



G[A] 



1 

7/7 + m 



777 — m 
(777) 2 — m 2 



(m-'yn)ijdse- i ^ m2 -^ n)2 l (2.7) 



0 

Obviously, the proper-time representation simply expresses an inverse opera- 
tor in terms of a convenient exponential function. Convergence of the proper- 
time integral at infinity is ensured by the implicit prescription m 2 — > m 2 — ie. 

1 Spinor indices will be suppressed in the calculation. 
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We are interested in the proper-time representation of the effective action 
W G K Hence, we have to solve the equation (cf. (2.2) and (2.4)) 

= - etr 7b M G ( x - X \ A )\ ( 2 - 8 ) 

for We now show that the ansatz 



iW (1) = i f d 4 x£ (1) = -i 



d s 
s 



tr x 



e is(7J7) 2 



(2.9) 



0 



fulfills (2.8) and thereby gives to within a constant. This constant has to 
be chosen so that the action satisfies the boundary condition that it vanishes 
for vanishing external field. For the calculation of the functional derivative of 
(2.9), we note that SII a (x ) = —eSA a (x): 



.SW^[A) 
1 SA a (x) 



oo 

I 7 jy-'-V fd‘yis VT, S -E!MEM .n-n- 



ie / dse _lsm tr 



sn a (x) 
7 a {x\^ne lsGn ^\x) 



^etr 



OO 

7 a {x\(m — 7 i7) i J d se - is[m2 - (7/T)2 



— etr[7 a G(x, x|H)], 



where in the third line we have made use of the fact that the trace of an odd 
number of 7’s vanishes. 

In this way, we have demonstrated the validity of (2.9), and so we can 
write for the unrenormalized Lagrangian 

OO 

£ (1) = [ — e - ism2 (x\e isGn)2 \x). (2.10) 

Z J s 
0 

For the propagator G(x,x'\A), as well as for the effective Lagrangian, we 
need to evaluate the object 

K(x, x'- s\A) = (x\ e isGn)2 \x'). (2.11) 



For the remainder of this section, we shall mainly be concerned with the 
various methods of finding an explicit representation in configuration or mo- 
mentum space. 

The first method that we shall discuss was proposed by Schwinger. It is 
based on the interpretation of K(x,x'\ s|H) as the coordinate representation 
of the proper-time “evolution operator” 
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U(s) = e~ iHs , (2.12) 

where we have introduced the “Hamiltonian” 

H = -( 7 JT) 2 = n 2 - n it = pp-eA lt , (2.13) 

and (J [xu = i/2 [7^,71/]. In this sense, the desired object K(x,x'\ s|A) is equal 
to a transformation amplitude (x| U(s) \x ') = {x, slx',0), while the evolution 
operator satisfies a “Schrodinger equation” 

i d s U(s) = HU(s). (2.14) 

To obtain the transformation amplitude, we have to solve the dynamical 
problem formulated by the following one-particle equations of (proper-time) 
motion of Heisenberg type: 

^ if] =2i7 M , 

( ^f = -i [i7„, H\ = e (F llv n v + n v F av ) + ( -a Xv d^F Xv 

= 2e F llv n v - ie d v F^ + ^a Xv d^F Xv . (2.15) 

Our strategy is as follows: first, solve the (operator) equations of motion 
(2.15) for n(s) as a function of x(s) and x(0); secondly, insert these findings 
into the Hamiltonian; and thirdly, solve the Schrodinger equation (2.14) in 
the coordinate representation, in which the Hamiltonian is now diagonal (for 
a pedagogical review, see e.g. [106]). 

Of course, this procedure can only be performed analytically for a cou- 
ple of simple external-field configurations, e.g. for constant and plane wave 
fields [148]. In Sect. 2.3, a detailed example of such a type of calculation is 
presented. 

The second method to be discussed is based on the classical WKB ap- 
proximation (see e.g. [147]). Here, we want to compute the transformation 
amplitude by means of 

(x', six", 0) = (27ri)- d/2 VDe lS , (2.16) 

where d = 4 is the dimension and h = 1, and the Van Vleck determinant is 

= 4 <- 1)4 det ; (2d7) 

S denotes the classical action. For the constant-field case, this was treated in 
[57], and also in [97] and [ 55] using the heat kernel method in QED (with not 
necessarily constant fields). As a side remark, we note that the transformation 
amplitude can also be expressed by means of the path integral representation: 

)= J Vx(X)e iS ^ x) \ 

x(0 )=x" ,x(s)=x' 



(x' , s|x", 0) = K{x ' , x"; s 



(2.18) 
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where 

S 

S = J dXL[x(X),x(X)}. (2.19) 

0 

The Lagrangian associated with the Hamiltonian of (2.13) is given by 

L = (2.20) 

Now, we have to calculate the classical action (2.19). This can be achieved by 
first solving the equations of motion (2.15), which, for constant fields, reads 
(77 = x/2 = p — eA) 

x{X) = 2eFi(A), (2.21) 

where we have used matrix notation, i.e. Fx = F ti v x v . This equation can be 
solved with the ansatz 



a;(A) = e 2eFA i(0), (2.22) 

and a further integration, together with the initial conditions x(X = 0) = x" , 
x(X = s) = x', yields 

1 p 2 eFA 1 

i(0) = ^ FirT 2eF(x'-x"), x(X) - x(0) = -^^(x'-x"). (2.23) 

This information is sufficient to compute the classical action, which turns out 
to be 

/ 

X 

S(x’,x";s) = e J d^A»(0 

x" 

+ -j(a/ — x") a eF a 13 [coth(eFs)] ^[x' — x") 7 + ^a^F^s. 

The last step on the way to evaluating (2.16) involves two partial derivatives 
of S with respect to x' and x" (cf. (2.17)), which at last brings us to the 
desired quantity, 

/— i ( 1 sinheFsA 

' /S= S5 exp ( - 2 trln ^r-)- 

According to (2.16), the complete transformation amplitude for constant 
fields is given by the result 



167T 2 s 2 
x exp 



1 , sinheFsA 



(x\ s\x" , 0) = —7^35 3 exp ( — ^tr In : ^ 



- ( 2 / — x") ( e F coth e Fs) (x' — x") 



) exp \ ie Jd£ tl A IJ '(g)j 

e i(e/2) CT F S) (2.24) 
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where the entire gauge dependence has been isolated in the form of the holon- 

omy factor (counter-gauge factor) exp^ie J d^A^)^, and the integration 

is along the straight path connecting x' and x" . 

A useful result is given by (cf. (2.18)) 



J Vx(X) exp 

x(0)=x" , x(s)—x' 



dA 



1 

4 X ^ 



ie / ^ / 1 sinh eFs 

= e ®" n — exp — tr In 



167T 2 S 2 



eFs 



x exp 



Ax' — a/^eFcoth eFs)(a/ — x") 



(2.25) 



Finally, we shall present a third method to evaluate the transformation 
amplitude K(x,x'\ s|A) that, on the one hand, represents a more pedestrian 
way of solving the problem but, on the other hand, supplies us with the 
Fourier transform of K{x, x ' ; s|A). The latter turns out to be extremely useful 
for finite-temperature calculations. 

In the following, we confine ourselves again to constant-field configurations 
where F =const. and work in the Schwinger-Fock gauge 

A ^ = - l -F^{x-x') v . (2.26) 

First, we move one step backwards and consider the object 



A(x, x'\A) = i / dse im s K(x,x'; s\A), 



(2.27) 



which is contained in the propagator of a Dirac particle (in coordinate rep- 
resentation) (2.7). From the Green’s function equation (2.5) for the Dirac 
particle, we can easily read off the corresponding Green’s function equation 
for A(x, aj'l-Asr) in the Schwinger-Fock gauge: 



S(x 



x') 



(2 = 5) 



( 7 U + m) G(x, aj'l-Asp) (2 = 7) [m 2 - (qTI) 2 ] A(x, oj'IAsf) 
m 2 — d 2 — e -{x-x')^F^F” x )(x-x')„ (2.28) 






^(xjx'IAsf)- 



Translational invariance reduces the coordinate dependence of A: A = 
A(x — x'|Asf)- The differential equation (2.28) is further simplified by ob- 
serving that any term except for the last one in the equation is invariant 
under (generalized) Lorentz rotations. Since the differential operator acting 
on A(x — 2 :'|Asf) contains the generator of these rotations in this last term, 
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namely — x v d^), A(x — ic'I^sf) must also be invariant. We conclude 

that 



-F^(x^d v - x v d„) A(x - x'I^sf) = 0. (2.29) 

This leads us to the problem of solving the Green’s function equation 



- d 2 + m 2 - - e -( x - x ')^ x F" x (x-x') v 



A{x - x'I^sf) 

= 8{x—x'). (2.30) 



It is helpful to study this equation in Fourier space by introducing the mo- 
mentum description 

A(x - x'\A SF ) = j -0^ e 1 **-*') ^(p|A sf ). (2.31) 

The differential equation (2.30) can then be written as 

( p2 + k2 + z1(!,| ' 4sf) = (2 ' 32) 

where we have used the short form k 2 = m 2 — ( e/2)aF . 

The solution is given by [55] 

oo ioo 

Z\(p|A SF ) =i J dse- iK - 2s e~ M{is) = J dse- K2s e~ Mis \ (2.33) 

o o 

where, in the last step, we have performed the substitution is — > s (no contour 
rotation!), and M(s) can be decomposed according to 

M(s) = p a X a ^(a)p p + Y(s). (2.34) 

The quantities X and Y depend additionally on the field strength and are 
given by 

Y(s) = -tr ln[cos(eFs)], (2.35) 

w , ^ tan(eFs) , 0 oc , 

X(s) = — , (2.36) 

where we have used matrix notation, e.g. F = (F)^". 

From the definition of A(x — x'|Asf) ((2-27) and (2.31)), we can read off 
the expression for the transformation amplitude K(x — x': s|71sf): 

K(X-X'-S\A SF ) = [ ^-P e ^-Y) e i(e/2)aFs e -Y(is) e -pX(is)p_ ( 2 .37) 

J (27r) 4 

For the purpose of obtaining the Lagrangian, we are especially interested in 
the coincidence limit x' — > x. We obtain 
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K(s\A SF ) = K(0-s\A SF ) 



d P i(e/2)aFs -Y (is) -pX(is)p 



{2fY 



(47t) ; 



e i(e/2) CTj F s e -Y( is) [ det X(is)]- 1/2 



(2 = 5) * e Ke/2),Fs 

(47t) 2 



— -tr lncosh(eFs) 



1 /IX 

x 7 r N2 ex P “o tr ln — 
(is) 2 \ 2 is 



(2.36) 



(47rs) 2 
= (xle 



1(7 n) s |x), 



(2.38) 



•'( exp ( — hr ln sillll(eFs) ) 
2 eFs ) 



(2.39) 



which is the result we obtained in (2.24) for x' = x" = x. To complete 
our study of the proper-time transition amplitude, the question of gauge 
dependence has to be discussed. Again, we start with the Green’s function 
equation (2.5) in coordinate space for an arbitrary choice of gauge: 



7 (— id — eA) + to 



G(x, x'\ A) = 5{x 



x'). 



(2.40) 



Replacing A M by a different gauge A ' = A^ — d^A (where A(x) is continuously 
differentiable), this equation reads 



_ e -i eA(x) 



— g-i eA(x') 



5(x — x') = y(— i d—eA + edA) + m G{x,x'\A — dA) 

7 (-id - eA) + ml e ieA{x) G(x, x’\ A - dA) 

7 (-id - eA) + m]e ieA(x) G(x, x'\ A- dA) 

7 (-id - eA) + ml e ie [ A W-MA)\ G{x,rf\A- d A), (2.41) 



where, in the third line, we have made use of the properties of the S func- 
tion. Comparing (2.41) with (2.40), we can identify the gauge transformation 
property of the Green’s function: 

G(x,x'\A-dA) = g-i^G)-^')] G {x,x'\A). (2.42) 

Since A(x) can be viewed as a “potential” for the difference between two 
equivalent choices of gauge, we may write for (2.42) 

iefd^A'^O-A^o] 

G(x,x'\A)=e x ' G(x,x \A). (2.43) 

Since A{x,x’\A) and K(x,x'-,s\A) are linearly related to G{x,x’\A), they 
behave equivalently under gauge transformations. 
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Taking these considerations into account, the proper-time transition am- 
plitude in an arbitrary gauge can be obtained from that in the Schwinger- 
Fock gauge according to 



K{x , x 7 ; s|A) = F{x, x'\ A) K(x , x'\ s|Asf), 
where we have introduced 

<P(x, x'\A) = exp 



X 

ie J df, (A*® + x')„) 

x' 



(2.44) 



(2.45) 



Since <£(x,x' | A) represents a mapping from a path in coordinate space onto 
the gauge group U{ 1), it is called the holonomy factor (or simply the holon- 
omy). Note that ^(x, x'| A) is invariant under continuous deformations of the 
integration path, because the integrand is curl-free. 2 In the coincidence limit 
x' — ► x, the holonomy factor simply reduces to 1, which explicitly demon- 
strates the gauge invariance of the effective action. 

This completes our search for an appropriate expression for the transfor- 
mation amplitude K{x — x'\s\A) or its coincidence limit K(s\A). In (2.39), 
we have rediscovered the findings of Schwinger [148]. The results of (2.37) 
and (2.39) are sufficient to finally evaluate the propagator of a Dirac particle 
and the effective Lagrangian according to the formulas (2.7) and (2.10). 

We conclude this section by stating the final expressions for the desired 
objects. The Green’s function for a Dirac particle (the electron propagator) 
in a constant external electromagnetic field (in an arbitrary gauge) reads 



G(x,x' | A) = $(x,x'\A) 



OO 

1 f ds 



(47t) 2 



m - + eF] M1/ (x-x') 1 



xe _ im s—L(s) exp 



-(x— x , )f(s)(x— x') 



i(e/2 )aF s 



where we have used the abbreviations 



r/ . ._ . . _ . , . . 1 sinh(eF.s) 

f(s) = eF coth(eFs) and L(s) = -tr in — . 



(2.46) 



(2.47) 



and 



ie/d«4A'‘(5)+(l/2)F^(5-x')d 

^(x, x|A)=e x ' (2.48) 

carries completely the gauge dependence of the propagator. 

According to (2.10) and (2.39), the one-loop effective Lagrangian for con- 
stant external fields can be written as 



These statements are not strictly tied to the use of the Schwinger-Fock gauge. 



2 
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OO 

1 f ds 



(es) 2 \G \ cot es{\/T 2 + G 2 +T) 1/2 (2.49) 



x coth 



esiy/F 2 + g 2 -T) 1/2 +-(es) 2 JE_l 

o 



where we have introduced the gauge and Lorentz invariants of the electro- 
magnetic field, 



T = (2.50) 

g = *F **" = ^ vX °F Xa . (2.51) 



The last two terms in the curly brackets of (2.49) stem from a field strength 
and charge renormalization. 

An expansion of the integral for small values of e corresponds to a weak- 
field approximation F^ v -C m 2 and yields 






° u T7 2 I _ 

45 m 4 45 m 4 

[(.E 2 - J3 2 ) 2 + 7(E ■ B ) 2 ] . 



(2.52) 



Euler, Kockel and Heisenberg [75, 98] and, independently, Weisskopf [ 
were the first to discuss effective Lagrangians of this type. 

Equation (2.52) can be interpreted as the first nonlinear correction term 
to the Maxwell Lagrangian, which defines the linear theory of classical elec- 
trodynamics. 

This concludes our introductory remarks about the proper-time method, 
which turns out to be a powerful tool with which to analyze the underlying 
(fundamental) gauge field theory of quantum electrodynamics. 



2.2 Polarization Tensor in External Fields 

We derive a new representation of the one-loop photon polarization tensor cou- 
pling to all orders to an arbitrary constant electromagnetic field, fully maintaining 
the dependence on the complete set of invariants. The problem of extensive Dirac 
algebraic computations is avoided by constructing a one-to-one mapping from the 
special case of parallel electric and magnetic fields, as treated by Urrutia, to the gen- 
eral case. The polarization tensor is decomposed into its eigenvalues and eigenspace 
projectors. We work out explicitly certain limiting cases, e.g. purely magnetic fields, 
crossed fields and arbitrary weak electromagnetic fields. Finally, the modification 
of light propagation in external fields is discussed. 

The polarization tensor 77 /iy (fc), also known as the photon self-energy, can 
be viewed as a modification of the (classically trivial) photon structure by 
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virtual processes of vacuum polarization. To one-loop order, it is an electron- 
positron pair which comes into existence for a certain time and annihilates 
itself during the propagation of the photon. 

On a heuristic level, the physical consequences of these vacuum polariza- 
tion effects can be illustrated by an intuitive picture: the virtual transitions, 
e.g. the electron-positron loop, confer the properties of the participant par- 
ticles on the photon. By this means, the photon can therefore acquire, for 
instance, a “size” (of the order of the Compton wavelength of the participant 
particle) or a “charge distribution” . 

In the trivial vacuum, i.e. in the absence of external perturbations, these 
internal photon properties are hidden; the photon propagates with the speed 
of light on the trivial light cone , which is defined by k 2 = — (k 0 ) 2 + ( k ) 2 = 0, 
where denotes the wave vector of a photon of arbitrary frequency. 

This situation changes substantially in the presence of various external 
modifications of the vacuum, even if an interaction with the propagating light 
is classically forbidden. These classically forbidden processes are, of course, 
the most interesting ones, and will be elucidated in this and the following 
sections. Within our heuristic picture, the vacuum modifications can interact 
with those properties of the photon which are caused by the virtual particles. 

In the present section, we consider constant, but otherwise arbitrary, elec- 
tromagnetic fields as the external perturbation. The study of the polarization 
tensor II fJ ’ u (k\A) will therefore serve as an investigation of the properties of 
the QED vacuum with external fields. The derivation of this kind of polariza- 
tion tensor has a rather long history. The properties of the tensor in a constant 
magnetic field have been studied by many authors [14, 124, 45, 159, 47, 152]; 
focusing on light propagation, a comprehensive investigation was performed 
by Tsai and Erber [160]. Another special case, the crossed-fields configuration, 
has been investigated by Narozhnyi [126] and Ritus [139]. However, a gen- 
eralization to arbitrary constant electromagnetic fields in a straightforward 
manner is associated with a substantial increase in calculational difficulties; 
detailed results like those for the purely magnetic case have still not been 
worked out. Although Batalin and Shabad had already obtained a represen- 
tation in proper-time form for the polarization tensor in an electromagnetic 
background in 1971 [22], their extensive result was only brought into a practi- 
cal form almost 20 years later [11]. A calculation employing special operator 
techniques has been performed by Baier et al. [15]. Within the world-line 
formalism, the polarization tensor has recently been calculated by Schubert 
[146] with comparably little computational effort. 

In the following, we shall derive the polarization tensor in an arbitrary 
constant electromagnetic field in a simplified manner; our approach is based 
on the findings of Urrutia [162], who solved the problem for the special case of 
parallel electric and magnetic fields. We shall demonstrate that the general 
case can be deduced from this particular one by an appropriate Lorentz 
transformation. Applications will be described afterwards. 
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2.2.1 Derivation of IT P" 



To begin with, we shall outline the direct calculation of 77 A11/ and point out 
where one encounters technical difficulties. According to the usual rules of 
perturbation theory, the polarization tensor can be calculated from 

n^(k\A) = — ie 2 tr 7 J 7m dip) 7^ gip ~ k), (2.53) 

where we have chosen to work in momentum space. Therefore, g(j>) denotes 
the Fourier transform of the Dirac fermion propagator, which, according to 
(2.7) and (2.27), can be written generally as 

G(x, x'\ A 1 ) = x'\ A' — A) (m — ift) A(x, x'\ A). (2-54) 



In the following, we shall use exclusively the Schwinger-Fock gauge for con- 
stant fields, Ag F = — (l/2)F liV (x — x') v , which allows us to express the gauge 
potential in terms of the field strength as well as to forget about the holon- 
omy factor <j)(x, x'\A!—A). Introducing the Fourier representation Z\(p|Asf) of 
A(x, x'|Asf), we can easily derive the Fourier representation of G(x , z'IAsf): 

d A p 



G{x, x / |A S f) = 



( 2^) 4 ~ 

[ d4p 

I (2tt) 4 



a i p(x-x') 



m-J a [p + 



l lp aK — 

2 dp K 



a i p(x-x') 



g{p)- 



Z\(p|A SF ) 

(2.55) 



In (2.33)-(2.36), we obtained the explicit form of Z\(p|Asf); inserting these 
results into g(p) and performing the momentum differentiation leads us to 



g(p) = i j d se-” 
o 



s p- Y ( is ) 



m — 7 a (p — ieFXp) c 



e -pXp e i(e/2)so-F 



(2.56) 



where we have employed matrix notation, and the quantities Y and X are de- 
fined in (2.35) and (2.36). After insertion of (2.56) into (2.53), the remaining 
problems are some Gaussian integrations and various Dirac traces; in fact, the 
Dirac traces can be looked up in Appendix C ((C.25)-(C.31)); these obviously 
demonstrate that the direct computation will be more than tedious. But the 
main problem is the reorganization of the terms at the end of the calcula- 
tion, especially the simplification of the trigonometric and hyperbolic func- 
tions. The latter would even prohibit successful use of algebraic-manipulation 
programs, at least for higher-order calculations. 

As already mentioned, the calculation becomes enormously simplified for 
a selected class of field configurations, such as the case of parallel magnetic 
and electric fields. Since space-time can effectively be decomposed into a 
transverse and a longitudinal subspace with respect to the field axis, the 
four-dimensional tensor analysis reduces to two-dimensional. If we choose the 
field to point along the 0 axis, without loss of generality, a vector decomposes 
according to 
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W + k^, k^ = (fc°,0,0,7 3 ), jfe£ = ( 0,fc\fc 2 ,0). (2.57) 

In the same manner, tensors can be decomposed, e.g. g . With 
respect to each subspace, we easily find the unique vector orthogonal to a 
given vector: 

fcfj* = (fc 3 ,0,0,fc°), k^ = (0,k 2 ,-k 1 ,0). (2.58) 

The polarization tensor for this special field configuration has been calculated 
in [162] in the way outlined above. We simply state the result here: 



0 




zz' 

sin z sinh z' 



(g^k 2 - k^k v )No + {g^k 2 - k^N, 

+ (<?f k i ~ K k ±)N 2 - + kft v ± )N 3 



(2.59) 
+ c.t. > , 



where “c.t.” means “contact term”. The electric and magnetic field strengths 
E, B are contained in the variables z = eBs and z' = eEs. The exponent <j > q 
is given by 3 



(j>o = m 2 + 



k \\ cosh z' — cosh vz' 
2 z' sinh z' 



The functions TV,: are 



k\ cos vz — cos z 
2 z sin x 



(2.60) 



No = cosh vz cos vz — sinh vz sin vz cot z coth z , 

„ t „ cosh z' — cosh vz' 

Ni = 2 cos 2 — — N 0 = TVi - N 0 , 

sinh z 

„ r „ , , cos vz — cos z 

N 2 = 2 cosh z N 0 = N 2 — N 0 , 



sin 2 z 



(2.61) 



N* = 



1 — cos z cos vz cosh vz' cosh z' — 1 



sin Z 



sinh z' 



— sin vz sinh vz' , 



where we have also defined the functions 7 Vi, 2 for later use. The contact term 
is, finally, determined by the requirement that 77 should vanish in the limit 
of vanishing fields and for wave vectors on the mass shell (light cone), k 2 = 0. 
This corresponds to a charge and field strength renormalization and ensures 
that the real photon will remain massless. We obtain 

c.t. = - e - im2s (l - v 2 )(g' ll 'k 2 - k^k v ). (2.62) 



Now, our claim is that the Lorentz-invariant form of the polarization ten- 
sor for an arbitrary constant electromagnetic field can be completely recon- 
structed from the special form given above for parallel electric and magnetic 



This formula has been misprinted in [162]. 
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fields. The situation should be familiar from the simpler case of a problem 
where it is known that the solution depends on the field invariants only. In 
that case it is sufficient to solve the problem for the special case of parallel E 
and B fields, i.e. to go to a Lorentz frame in which E\\ ± B, and reconstruct 
the invariants later on according to T = (1/2) (B 2 — 77 2 ) and Q = — EB . This 
procedure leads to the correct results, since there is a one-to-one correspon- 
dence between the two sets of linear independent variables 77, B — > J-,Q. 

Returning to the present problem, we have to establish such a one-to-one 
correspondence between the dynamical variables contained in (2.59) and a 
complete set of linearly independent Lorentz and gauge invariants, as well 
as a one-to-one correspondence between the tensor structure of (2.59) and a 
complete tensor basis for the polarization tensor in arbitrary fields. 

Let us begin with the scalar quantities. In the present problem, the basic 
building blocks for the construction of gauge- and Lorentz-invariant scalars 
or pseudoscalars are given by the field strength tensors F ^ and *F tJl ' y and 
the wave vector k The number of linearly independent invariants is four, 
as can be proved by a simple group-theoretical argument (cf. Sect. 3.5). Here 
it is convenient to choose the secular invariants a, b as defined in Appendix 
B for the purely field-dependent invariants, since they reduce to a — > B and 
b — > —E in the case of antiparallel electric and magnetic fields. Note that 
we shall assume E < 0 in the following without loss of generality to obtain 
a simple invariant structure (cf. Appendix B). A convenient choice for the 
remaining invariants is given by 

Zk=(k a F aK )(k p F^ K ), (2.63) 

1,2 _ jLAtl, 

lb lb IV f_L • 

In fact, there are exactly four dynamical variables which 77 in (2.59) de- 
pends on: 7?, B, k 2 and k\. The map from these variables onto the set of 
invariants is simply given by 



a = 77, b = —77, 

z k = - E 2 k\ + B 2 k 2 ± , k 2 = fcfj + k 2 ± , (2.64) 

which is valid in a system where B\\ — E. The inverse map is obtained by a 
simple calculation; the nontrivial relations are 



2 _ a 2 k 2 - z k 
II a 2 + b 2 



k 2 ± = 



b 2 k 2 + Zk 
a 2 + b 2 



(2.65) 



Hence, the one-to-one correspondence between the scalar quantities is estab- 
lished. 

Let us turn to the tensor structure of 77 ^ . There are exactly four lin- 
early independent vectors which we can construct from and F we shall 
employ the choice 

W, Fk ^ = F^ a k a , FV = F^F^k 0 , *Fk» = *F^ a k a . (2.66) 
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(Note that vectors containing higher powers of the field strength tensors can 
be decomposed with the aid of the fundamental algebraic relations (B.3), 
(B.4).) Since 77 is a symmetric second-rank tensor by construction, it can 
only possess ten independent components; indeed, we can only construct ten 
independent symmetric Lorentz second-rank tensors from the four vectors 
given above. 

Next, we take advantage of Furry’s theorem, which in this case tells us 
that the polarization tensor can only contain an even number of powers of 
F in each term. This is because, at each external line, a field strength tensor 
is attached to the loop, and there is only an even number of external lines. 
Since the invariants, as chosen above for example, only involve even numbers 
of F ^ (note that a and b are defined via and the possible 

Lorentz tensors must also contain only even numbers of F^ . This leads to 
the precept that we are only allowed to combine the four vectors in (2.66) 
in such a way that the tensor product contains even powers of F^ . So, the 
number of possible tensors reduces to six. 

Finally, we employ the Ward identity as a consequence of gauge symmetry, 
which demands that n^k u = 0. This reduces the number of basic tensors 
to four. In other words, there are only four linearly independent tensors that 
generate the tensor structure of Ft fJ,l '(k\A). In fact, the polarization tensor for 
the special case (2.59) is constructed from four linearly independent tensor 
quantities which respect all of these symmetries: 

g^k 2 -Wk\ g^q-k^kl g^kl-k^kl, k^ + kfk v ± .( 2.67) 

Our final task is to establish a one-to-one correspondence between these ten- 
sors and an appropriate set of Lorentz-invariant tensors constructed from the 
vectors (2.66) in the special Lorentz frame where B | — E. 

As a first step, it is useful to note the following identities: 

g ii *u *n *ii - *n *ih 

g^ki-k^kl = k^kl, ( 2 . 68 ) 

which can be verified by a straightforward calculation employing the defi- 
nitions (2.57) and (2.58). Next, we represent the Lorentz-invariant vectors 
(2.66) by and k^ ± in the special Lorentz frame. By simple matrix mul- 
tiplication, we find 

Fk M = Ek^ + Bk^ = -bkX + ak^, 

*Fk fJ ‘ = Bk fj* - Eh%_ = afcfj* + bk 

F\^ = E 2 fcff -B 2 k^ = b 2 fcff - a 2 k^, 

r ./ 4 — r-A 4 _i_ r-A 4 

iXj — ruy | . 

These equations can easily be inverted: 



(2.69) 
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^ 4 «=') • = p^p ( rt ' + “ 2 ' 

•Z = p^p (* + « ra ") . = p^p (- rt “ + !> 2 *") ■ (2-70) 



which are valid in the special system. However, the right-hand sides of (2.70) 
are written in a Lorentz-invariant form and can therefore be transformed into 
any Lorentz system. For instance, let us define the Lorentz transformations 
of kt? ± as 



~ k \\ - 


^11 a 2 + b 2 1 


(a*Fk» -bFk 11 ) 


K ) 




1 




- 


^ a 2 + b 2 1 


[b*Fk» + aFkn 



(2.71) 



where the subscripts || and _L are intended to remind us of the meaning of uy 
and u_l in the special Lorentz system (the longitudinal and transverse parts 
of fc). As a cross-check, we can prove the following identities in a general 
system; they are, of course, also valid for kJ^ ± in the special system: 

v l = v ii v t>‘ = - a ^+p k ( s *if = _t i)- 

(= ft! = fci). (2-72) 

Wy 1 V±tJ, = 0 (= k Ij* k±^ . 

Inserting the Lorentzian generalizations of k^ ± and fcy 1 ± given in (2.70) into 
the tensor quantities of the special system (2.67), we find a Lorentz-invariant 
tensor basis for the polarization tensor in an arbitrary constant electromag- 
netic field: 



n^k 2 - k^k 1 ' 

9 ii K \\ nrii 
sffci - k^kl 
kY~k\\ + ~kfil 



U 1/ 

«j_ W -L> 

U jy - U 1/ 

v l v \\ +< Pu 



g^k 2 - k»k v -> k 2 (Pg u + Pjj*" + P £") , 



where the projectors P ^ jj ± are defined by 



(2.73) 



pr = 



k 2 



2Pz k /k 2 +Q 2 -{z k /k 2 f 






P^ v = 



P^ = 



i>7 v 



_L ^_L 



(2.74) 



and satisfy the usual projector identities Pq y j_ = Po,||.±, -Po,U,i% = 1. 
Moreover, the projectors are orthogonal to each other: Pk 1 ' Pj = 0 for 
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i 7 ^ j and i,j = 0, ||,_L. We may furthermore introduce the short form 
Q = v^v 1 ^ + which is neither a projector nor orthogonal to the 

iffl’s but is orthogonal to Pq 1 ' . 

We are finally in a position to transform the polarization tensor for the 
parallel field configuration into its generalized form for arbitrary constant 
electromagnetic fields : 

n^(k\A) = n 0 p£ v + 77 || p» v + n± p^ + oq^ v , ( 2 . 75 ) 

where i7 0 11 j_ and 0 are functions of the invariants and are given by 

oc 1 

Tin I a f ds f dv 

nlrSTjl 

{0} 01 

(2.76) 



/ 


0 


\ 


c -i -cko eaS ebs 


N 0 v 2 ± - Niv^ 


L, 


sin eas sinh ebs 


N 2 v 2 ± - N 0 vj 


V 


{ -n 3 J 


/ 



Substituting the invariants into (2.60) and (2.61), the functions jV, and 0o 
yield 



0 o = m 2 - 



' u ll cosh ebs — cosh ize&s cos veas — cos eas 



2 ebs sinh ebs 2 eas sin eas 

No = cosh vebs cos veas — sinh vebs sin veas cot eas coth ebs , 

~ cosh ebs — cosh vebs 

N i = 2 cos eas — 



N 2 = 2 cosh ebs 



sinh ebs 

cos z^eas — cos eas 

• 2 > 
sin eas 



(2.77) 



W = — 



1 — cos eas cos z'eas cosh z/efrs cosh ebs — 1 



sin eas 



sinh ebs 



- sin veas sinh vebs. 



The scalars ± are given by certain combinations of the invariants and can 
be found in (2.71). The contact term given in (2.62) contributes equally to 
the Ili s, 



c.t. = -e~ im2s k 2 (l - v 2 ), 



(2.78) 



but does not modify the function 0 , which is already finite. 

Although (2.75) is the most convenient representation, a physical meaning 
can only be ascribed to the particular terms in a diagonalized version. While 
Pq V indeed projects onto an eigenspace of IP 11 ' with eigenvalue 77o, this 
is generally not the case for the projectors owing to the mixed term 

0 . However, the function 0 vanishes in those cases where either a or b 

vanishes, i.e. E • B = 0, since A 3 (a = 0 or b = 0) = 0. Then, the vectors 
v 1 ^ ± indeed represent eigenvectors of with eigenspace projectors P^ ± 
and eigenvalues in other cases, the eigenvectors are formed by linear 

combinations of v\f and v ((. These can be found by a direct computation, 
which we shall encounter in the next subsections. 
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The fourth eigenvalue of II ^ is, of course, zero, since the remaining 
eigenspace projector P£ v = g ^ — Ptf 1 ' — P^ v — P = k^k v /k 2 satisfies 
the Ward identity i7 M , y (fc) P£ x = 0. 

In view of the technical problems associated with a general, direct com- 
putation of II ^ for arbitrary constant field configurations, the simple ten- 
sor formalism developed so far reduces the calculational efforts considerably. 
Moreover, the use of the one-to-one correspondences between the invariants 
and the basic tensors established here is not limited to the one-loop level, 
but remains applicable to all orders in perturbation theory. 

2.2.2 Applications 

This representation of the polarization tensor in an arbitrary constant elec- 
tromagnetic field can be reduced to various well-known representations for 
the limiting cases of certain field configurations. As a demonstration of the 
physical information contained in 77 A11/ , we shall apply our findings to the 
question of how the propagation of light is modified in an electromagnetic 
background. Here, we mainly focus on low-frequency photons in order to 
make contact with the effective-action approach in chapter 3. Nevertheless, 
the polarization tensor contains information about light propagation at any 
frequency scale. 

Having integrated out the fermions (electrons and positrons), the effects 
of vacuum polarization can be included in the effective action for the electro- 
magnetic field by writing 



Weff = f d 4 xC cS = - [ 

J 7 4 A+a 

— i J A i xd i x' a ll {x)F[ fJ ' u {x, x'\A)a v (x r ). (2.79) 



In view of the physical situation under consideration, we have separated the 
electromagnetic field into a constant background part A ^ and a plane wave 
field with potential a M (fc) and wave vector Deriving the equation of 
motion for the plane wave field in momentum space leads us to the light cone 
condition for the propagating wave: 

fcV" - W + n^(k\A)\ a v (k) = 0. (2.80) 

Employing our knowledge about the tensor structure of TP*", the light cone 
condition simplifies if the polarization of the plane wave lies entirely in an 
eigenspace of FI /i ". For example, if a M (fc) = P^a u (k), i.e. a At (fc) is propor- 
tional to F 2 k ^ + (zk/k 2 )^, then we obtain the light cone condition for a 
“longitudinal” photon, 

k 2 + n 0 = 0, 



for a M (fc) = PQ V a„{k). 



(2.81) 
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Since 77 0 ~ k 2 (cf. (2.76)), (2.81) can always be solved trivially by k 2 = 0. 
From this, we may infer that there are no longitudinal photons. The interest- 
ing question is whether there are further, nontrivial solutions to (2.81) which 
hint at the existence of longitudinal photons and a new mode of propagation. 
Different points of view on this subject can be found in [47], [152] and [161]. 

The existence of non-trivial solutions to the light cone condition should 
be expected if the polarization of the plane wave lies in the plane spanned 
by since these polarizations are continuously related to the transverse 
ones in the trivial vacuum. The light cone condition then reads 



(fc 2 + n n )pf + (k 2 + n ± )p£ v + e or 

for a M (fc) = {P? v + P^ l ')a u (k), 



,{k) = 0, 



(2.82) 



which we shall refer to as the light cone condition for transverse photons. Of 
course, the physical polarization modes still have to be identified by diago- 
nalization of (2.82). 

For the following applications, it is useful to be aware of some identities 
that relate the Lorentz-covariant objects to their Galilean counterparts in 
a certain reference frame; these identities can be checked by simple matrix 
algebra (cf. Appendix A): 

= (u>, k) = k(u, k), 

Fk v = k(k • E, k X B + vE), 

*Fk^ = k(fc • B, k x E + vB), (2.83) 

FV = k [E ■ (kxB) + vE 2 , (k • B)B + (fc • E)E - B 2 k + vExB\ 
Zk = k 2 [B 2 sin 2 6 b + E 2 sin 2 0 E - (1 - v 2 )E 2 + 2 vE • (k X B)] . 



In the first line, we have introduced the phase velocity of the plane wave, 
defined by v = tu/|fc|; we have furthermore employed the short form k = |fc| 
(the Lorentz scalar k 2 should not be confused with the 0(3) scalar k 2 = |fc| 2 ). 
In the last line, we have defined the angles 6b and d E between the magnetic 
and electric fields and the wave vector k via B 2 sin 2 6b = B 2 — (B • k) 2 and 
similarly for d E - In terms of the phase velocity, a light cone condition of the 
form k 2 + II = 0 can be written as 

v 2 = l + j|, (2.84) 

where 77 generally still depends on v 2 . 



Magnetic Field 

The problem of evaluating the polarization tensor becomes considerably sim- 
plified if at least one of the invariants vanishes. This is the case, for example, 
if a reference frame exists in which the electric field vanishes. We are left with 
a purely magnetic field, and the invariant structure can then be written as 



22 



2. Nonlinear Electrodynamics: Quantum Theory 



6 = 0, a = B , Zk = k 2 T> 2 sin 2 9b, k 2 = k 2 (1 — v 2 ) . (2.85) 



It is therefore useful to rewrite the problem in terms of 13, k, v instead of 
a,Zk,k 2 ; for convenience we shall omit the subscript B on the angle 9b- 
Furthermore the functions TV) and </>o reduce to 



Ng = cos veBs — v sin veBs cot eBs , 



N\ = (1 — i/ 2 ) cos eBs, 



N 2 = 2 



cos veBs — cos eBs 



( 2 . 86 ) 



4 >o = rn 2 — k 2 < [sin 2 9 — (1 — u 2 )] 



sin eBs 

1 — v 2 sin 2 9 cos veBs — cos eBs 
4 2 eBs sin eBs 



Since N 3 and, consequently, <9 vanish, 77||.j_ become eigenvalues of the polar- 
ization tensor, and the transverse light cone condition (2.82) decouples. For 
the different modes, we obtain 

k 2 + 77|| = 0, || mode : a M ~ 

k 2 + II± = 0, 1 mode : ~ v^_. (2.87) 

The eigenvectors for purely magnetic fields read: 

«jl =k [(k-B),vB], 

< = k(0 ,kxB). (2.88) 

Referring to (2.87), the vectors U|[ ± determine the polarization directions 
of the plane wave which diagonalize 77 A “'. The direction of the electric 
field vector of the propagating plane wave e is related to these vectors via 
ei = f° l = k°a l — k l a° ~ w vl ± — Wy ± k l . Hence, the _L mode is polarized per- 
pendicular to the plane spanned by the magnetic field B and the propagation 
direction k, while the || mode is polarized parallel to this plane: 

ey ~ v 2 B — (k ■ B) k, 

e± ~ k X B. (2.89) 



Note that the || mode is only “almost” transverse, owing to the factor of v 2 . 

Now we arrive at a subtle point with a long tradition in the literature: since 
we are dealing with a one-loop calculation, we expect the shift of the phase 
velocity to be small compared with the vacuum velocity. This assumption 
should at least hold for moderate magnetic fields 0 < B < 0(B CI = m 2 /e) 
and moderate frequencies 0 < u> < 0(m). In this case, we may set v = 1 
and k = u> in the polarization tensor, because a small deviation would lead 
to only a next-to-leading order correction to the velocity shift. Note that this 
assumption must be confirmed by consistent results. The first consequence is 
that the longitudinal light cone condition is already satisfied, since v = 1 <t=> 
k 2 = 0; in this regime, there is no longitudinal photon. The contact terms 
also vanish. The function <j > 0 reduces to 



2.2 Polarization Tensor in External Fields 



23 



(j> o = rrr — u> 2 sin 2 9 



1 — v 2 1 cos veBs — cos eBs\ 

4 2 eBs sin eBs ) 



(2.90) 



The scalars v jj j_ simply become v 2 ± = u> 2 sin 2 9. The relevant components 
of the polarization tensor can then be written as 



77 m j_ = xr sin 2 9 



27 r 



ds 



dv 



-is0o N 



,X> 



(2.91) 



where 



N» = 



eBs cos veBs 



N ± = - 



sin eBs 
eBs cos veBs 



— eBs cot eBs ( 1 — v 



sin veBs 



sin eBs J 



(2.92) 



veBs sin veBs cot eBs 



sin eBs sineTJs 

2e5s(cos veBs — cos eBs ) 
sin 3 eBs 



Here, we have arrived at the representation of the polarization tensor found 
by Tsai and Erber [160], who employed this expression to obtain refractive 
indices for a magnetic background over a wide range of frequency and field 
strength. We shall confine ourselves to the case of low-frequency photons, 
ui <C in, which allows us to approximate the function (f > q in the exponent by 
77o — > to 2 . Bearing in mind that we assume a small velocity shift, we may 
write for the phase velocity 



v = 1 — — sin 9 ) 
47T 



II, X, 



where Tj« x ar e functions only of the field strength: 

OO 1 



,_L — 



/ 



ds 

s 



dv 



— l m s AT 

' Wl.x- 



(2.93) 



(2.94) 



The necessary integrations can be performed analytically; details may be 
looked up in Appendix D (see (D.28) and (D.29)). We find 



B r 



77 



^ii ( fl ) = { ■ - 1 - of ) - 4 of ^ ( 0’ ) - ** ( i + of 



2 BJ 2 B 

/i cr B CI ( B C1 

+2^1n-^ -2)^ 
277 



B c 



Be 



2 B 



V±(B) = -4 
-2 



2 B 

Be r 
2 B 

B rr 



no 



Her 
2 B 



2 B 
4 



Her 
2 B 



2 B 

3 H cr 
2 

#11 + 



1 

3’ 

Her 

2H 



He 



In rki _ 2^1 - 4 
2H 2H 25 



Her 

25 



(2.95) 
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Our representation employs the first derivative of the Hurwitz zeta function 
('(s,q). With the aid of the identities of Appendix D (see (D.40)), its equiv- 
alence to the representation of [160], which involves a generalized gamma 
function, can be proved. 

In the limit of weak fields, B <C B CI , we recover the well-known result for 
the velocity shift: 



Tjil ~ 1 — ^—jB 2 sin 2 9 + 0 
" 45 to 4 

v_L ~ 1 - sin 2 9 + O 

45 m 4 



B 

B C i 

B 

B c , 



(2.96) 



The limit of strong fields can be studied by expanding (2.95) in terms of 
B CI /B <C 1. We obtain 



v\\ ~ 1 sin 2 9 

11 47T 

u_l — 1 — j- sin 2 9 
47T 



2 B 

3 B c , 



o 



1 -+ 2 1-,l 1 

3 3 



Act , B CI 
2 B ln 2 B 



O 



Bqt , B ct 

ln 

2 B 2 B 



(2.97) 



where 7 denotes Euler’s constant and L\ = 0.248754477 . . .. Note that the 
phase velocity for the _L mode remains close to one, even for strong magnetic 
fields. The maximum velocity shift amounts to duj_ max = — a/( 67 r) ~ —3.8 x 
10 -4 . This justifies the assumption regarding the smallness of the velocity 
shift at the beginning of the calculation for the _L mode. 

In contrast, the velocity shift for the || mode increases linearly with B 
for strong fields and seems to be unbounded. This, of course, contradicts the 
assumption that the velocity shift should be small; hence, (2.97) for the | 
mode is only valid as long as —5v± -C 1. This restricts the magnetic field 
strength so as to obey B/B CI < 7 r/a ~ 430. 

To obtain the correct strong- field behavior in the limit B — > 00 , one must 
evaluate the integrals without a restriction to small velocity shifts. In Sect. 3, 
we shall obtain the correct strong-field behaviour within the effective-action 
approach. In this approach, it is not necessary to perform tedious integrations. 

Let us finally remark that these phase velocities for the low-frequency 
limit do not explicitly depend on the frequency; hence the group velocity 
coincides with the phase velocity in this limit. It will turn out that this is a 
general feature of the low-frequency (= soft-photon) limit. 



Absorption in a Magnetic Field 

As an example of a high-frequency process, we investigate the absorption 
rate for a propagating photon in a magnetic field due to c + e _ pair produc- 
tion, following the work of Tsai and Erber [159]. By reference to the optical 
theorem, the absorption coefficients for the different polarization modes are 
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related to the imaginary part of the polarization tensor in a magnetic field 
(2.91) by 



ttii _l(w) = Im7T|| _i_. (2.98) 

Of course, it is hardly possible to evaluate the imaginary parts of ex- 
actly; therefore, we shall confine ourselves to the energy regimes relevant to 
the present situation. 

In particular, the exponential factor e _ll ^ oS has to be treated with great 
care. Substituting for the s integration by means of eBs = z, the exponent 
reads 



^° S = ~eB 



u“ . o „ t 1 cos uz — cos z 

z + z— sm 9 : 

^ '2 zsmz 



m 




(2.99) 



First, note that the term in the round brackets has the following expansion 
for small z: 



1 cos uz — cos z 

2 z sin z 




— (l-u 2 ) 2 z 2 
48 V 1 



0(z 4 ). 



(2.100) 



We are allowed to insert (2.100) into (2.99) if the principal contribution to 
the z integral stems from the lower bound, where z<l. Regarding (2.99), 
this is the case if m 2 / eB 1 and u/m sin# 1, i.e. for weak fields and high 

frequencies. 4 For photon absorption, the latter condition exactly matches our 
requirements, since pair production can only occur beyond the two-particle 
threshold, ui 2m. 

Therefore, we approximate the exponential factor by 

e -^° s ~e _iS , (2.101) 



where “ is defined as 



m 
= ~eB 



U 2 . 2 fl (l “ V 2 ) 2 3 

— - sm 2 6- — z A 



m * 
1 



48 



(2.102) 



In the last step, we have performed the substitution 



1 - v 2 



V = 



sin 9 



£ = 



8 m 2 m 1 



1 



3 eB lo sin 9 1 — u 2 



(2.103) 



For our further purposes, it is useful to introduce another parameter: 

4 1 



3 eB w 

A = - — — sm#, 
2 m z m 



and hence ^ = 



A 1- 



/ 2 ‘ 



(2.104) 



Throughout these considerations, it is understood that the z integral contour 
lies slightly below the positive real axis. 
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Note that the requirements m 2 /eB 1 and (w/m) sin0 3> 1 do not restrict 

A to a particular interval; A can become large for sufficiently high frequencies 
or small for weak magnetic fields. 

Since only small values of z are relevant for the integral, it suffices to 
approximate the functions N\\.± by the first terms of their Taylor expansions: 

N± = \( l-^Q + r 2 ) (2-105) 

Finally, the absorption coefficients can be written as 



OO 1 




where the representation of the modified Bessel function (Airy integral) 



1 

7! 



OO 

#2/3(0 = J dyy sin 










(2.107) 



has been employed. Equation (2.106) demonstrates that the absorption coeffi- 
cients are proportional to the Larmor frequency 10 l = eB /m. Apart from the 
angle dependence sin 0, the remaining formula is a pure number depending 
on the parameter A as defined in (2.104). 

In the various limits of A, one can make use of the usual asymptotic 
expansions of the Bessel function. For example, in the region where A > 1, 
i.e. sufficiently high frequency, we find, following [159], 



||, iM ^ \asmO— A" 1 / 3 ^F(3/2) 



x 



di/(l-i/ 2 ) 1/3 




0 



1 . 
= -asm 
2 




m 



1/3 



2 1 ' /3 v / 3 T ,2 (2/3) 
7^ #(7/6) 




[(3) || , (2)_l] 
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~ iasinfl— A 1/5 [(1.04)||, (0.69)j_] . 

2 TO 



(2.108) 



For small values of A -C 1, i.e. sufficiently weak magnetic fields, we obtain 






• n eB 
asm 0 — 

m 




X jdu( l-^-^e-I 4 /^ 1 - 2 ))! 



0 



1 • n eB 

-asm 0 — A 

2 mV 27rA 



OO [- 

x/d«(»-i)-^ (A + 5 1) 









_L 



e“ 4 “/\ 



(2.109) 



These findings coincide with the expressions found in [ 58, 117, 74, 14, 116, 
135, 166], 

Note that (2.108) and (2.109) describe the asymptotic absorptive behavior 
of the magnetic field. Both absorption coefficients tend to zero for A — > 0 or 
A — > oo. In between, a maximum exists at A ~ 16 for km/wl and at A ~ 19 
for k_l/wl. These features are exhibited in more detail in Fig. 2.1. 
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Fig. 2.1. Absorption coefficients kmj_ in units of the Larmor frequency cul = eB /m 
versus the parameter A as defined in (2.104); the propagation direction is assumed 
to be orthogonal to the magnetic field, i.e. 9 = ir/2, sin# = 1 



28 



2. Nonlinear Electrodynamics: Quantum Theory 



Crossed Fields 



A particularly interesting field configuration is provided by so-called “crossed” 
fields, where E and B are not only orthogonal but also of equal modu- 
lus: EEB and E 2 = B 2 . This causes the vanishing of the two invariants 
T = (1/2) (B 2 - E 2 ) = (1/2) (a 2 - b 2 ) and Q = E ■ B = ab. The field 
dependence is therefore completely carried by the invariant Zk = (k ll F fJ ’ u ) 2 . 

The structure of the polarization tensor in such a background was orig- 
inally studied by Narozhnyi [126] and Ritus [139]. In the present work, the 
crossed-field configuration serves as an instructive example of the question of 
how to take the correct limits for the invariants; it also represents a touch- 
stone for our formalism, since we derived our representation for II ^ from a 
completely different field configuration, which cannot be trivially related to 
the crossed-field situation by a smooth limiting process. 

Indeed, taking the limit a, b — > 0 requires special care. Problems arise 
from the fact that the choice a = b = 0 does not uniquely determine a class 
of Lorentz-equivalent field configurations; e.g. F = 0 implies a, b = 0, but 
a, b = 0 obviously does not force F ^ to vanish. Thus, we must ensure that 
our limiting procedure of taking a, b — > 0 does not automatically lead us to 
the trivial vacuum = 0. The latter would be the case if, for example, we 
first set b to zero and then took the limit a — > 0. Since 6 = 0 corresponds 
to the existence of a reference frame in which the field is purely magnetic, 
a = B, taking the limit a — > 0 would then definitely lead us to the trivial 
vacuum. The same undesired result would be achieved by interchanging the 
roles of a and 6. 

The correct procedure is to set a = b first, and then take the limit (a = 
6) — > 0. In terms of field strengths, we first choose \E\ = \B\ and afterwards 
rotate the direction of either E or B into the perpendicular position EEB. 
By construction, taking the limit (a = 6) — > 0 does not affect the electric 
or magnetic field strength but changes their orientation in a given reference 
frame. 

To obtain the polarization tensor in the limit of vanishing a and 6, it is 
sufficient to consider the functions </>o and N t in (2.77) for small values of a 
and 6: 



4>o 

N 0 

Ni 

N 2 

n 3 



6 



(es) 2 



(1 _ „2) + ^_(i _ i/ 2 )(— 36 2 + v 2 b 2 - 



(es) 2 



(1 — v 2 ) 2 ab. 



6a 2 ), 

6b 2 ), 



4 



(2.110) 
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Employing (2.72), we find the following identity in the limit (a = b) — > 0: 

N 2 v]_ - -/Voi| = k 2 ( 1 - v 2 ) + ^-(1 - iz 2 )(9 - v 2 ) z k 

= N 0 Vj_ — JViuji , for (a = 6) — > 0. (2.111) 

Referring to (2.76), we are able to conclude that the contributions to the 
polarization tensor 77|| and Il± are equal for crossed fields. Nevertheless, 
the polarization tensor is not degenerate, since the term 6 Q ^ does not 
vanish for crossed fields. 5 Therefore, the projectors Pjf 1 ^ are not related to 
the eigenspaces of W J ’ U . In order to diagonalize the polarization tensor, it is 
useful to note the following relation, which is valid for crossed fields (TTy = 
n± = 77 ||,_ l ) : 



n n pf + n ± p^ = n ]lx f^ + 



*Fk*Fk 

Zk 



ins 



( 2 . 112 ) 



where we have made use of the explicit form of the projectors P as de- 
fined in (2.71), (2.72) and (2.74). The contribution OQ 111 ' can similarly be 
represented; we find for the combination —N 3 Q^ U in (2.76) 



- N 3 Q^ = 



{esf 



(1 — v 2 ) 2 



Zk 



FkFk *Fk*Fk 



Zk 



Zk 



fiiy 



(2.113) 



In fact, the tensors Fk^Fk^/zk and *Fk fJ, *Fk v / Zk are orthonormal projec- 
tors in the case of crossed fields and they are orthogonal to the eigenspace 
projector Pq V . Hence, we arrive at the diagonalized form of the one-loop 
polarization tensor in crossed fields : 



= n 0 pr + Hi 



*Fk fJ, *Fk 1 ' 

Zk 



1 1> 



Fk^Fk" 



Zk 



(2.114) 



where 



Tin 



a f ds f dv 



n 1 =d- T( 1-Oe 



2 \ _ — im s 



n 2 



27 r 



k 2 + k fz k + 1 - u 2 )z k 

k 2 + ^fz k - te£(l - v 2 )z k , 

4>i = - v 2 ) k 2 + - v 2 Y z k - 



- k 2 



(2.115) 



For the derivation of (2.115), we have employed (2.111)-(2.113), which we 
have inserted into (2.76), taking the limit (a = b) — > 0. 

5 At first sight, the situation appears to be paradoxical: 0Q —> 0 for either 
a — * 0 or b — > 0, but QQ ^ remains finite for (a = b) — > 0. This demonstrates 
how important the correct treatment of the limiting process is. 
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The dynamical content of the polarization tensor (2.114) is completely 
determined by the two dynamical variables k 2 and Zk, while the orientation 
of the eigenspaces can be read off from (2.83) in the limit \E\ = \B\ and 

E ■ B = 0. 

A different representation of 77 ^ for the case of crossed fields, in terms 
of integrals over Airy functions, can be found in [139]; the corresponding for- 
mulas for the functions 77, are particularly useful for studying the analytical 
properties of 77^" in the complex k 2 plane. 

Finally, we apply these findings to the question of light propagation in a 
crossed-fields background; for convenience, we confine ourselves to the soft- 
photon limit ui <C m and consider weak fields only. Regarding Zk in (2.83), 
the weak-field condition can be summarized as the requirement for Zk/ k 2 to 
be smaller than m 4 : Zk/ k 2 <C to 4 . Again, we have to impose the additional 
assumption that the velocity shift remains small in order to set k 2 = 0 in the 
polarization tensor. This implies that we can neglect the function </i in the 
exponent of (2.115): e -18 ^ 1 — > 1. As an initial consequence, we observe that 
77o vanishes in this limit; this is not astonishing, since new longitudinal modes 
of photon propagation are not expected in the classical limit of weak fields. 
The remaining integrals in 77^2 can be carried out elementarily, leading to 
the following simple result to first order in Zk' 



7 ?! 

77 2 



= -2 



f 14/45 1 a 2 
{ 8/45 / m l 



Zk- 



(2.116) 



Upon insertion of this into the light cone condition (2.84), we find for the 
phase and group velocities of the respective polarization modes 



14 a 2 Zk 
45 to 4 k 2 ’ 
8 a 2 Zk 
45 to 4 k 2 



(2.117) 



Equations (2.117) are identical to the findings of Narozhnyi [126], who also 
studied the strong-field case. 

From the form of the eigenvectors Fk M and *T7e M as given in (2.83), we can 
conclude, similarly to (2.89), that the polarization directions of the electric 
field in the two modes are ( v 2 ~ 1) 

ei ~ — k X E + B — (B ■ k)B, 

e 2 ~ k X B + E - {E ■ k)E. (2.118) 



It is also useful to consider Zk / k in a given reference frame; for this, we also 
retain the assumption v ~ 1, and conclude from the last line of (2.83) that 

% = B 2 sin 2 0 B + E 2 sin 2 0 E + 2k • (B X E) 

= B 2 sin 2 Ob + E 2 sin 2 0 E — 2k • S 

= f 2 (l + (fc-n)) 2 . 



(2.119) 

(2.120) 
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Equation (2.119) expresses Zk/ k 2 in terms of the transverse electric and mag- 
netic fields and the Poynting vector S — E X B. Equation (2.120) explicitly 
exploits the properties of the crossed-fields configuration: we have defined 
a variable £ for the modulus of the field strengths £ = E = B and intro- 
duced the vector h normal to the plane spanned by B and E: h — B X E. 
From this representation, we can read off that a propagation direction in the 
crossed-fields case exists which renders no modifications: if k • h = — 1, then 
Zk = 0 = Sv i j2 . In other words, the propagation of light is not affected by 
weak crossed background fields if (k, E, B) form the basis of an orthonor- 
malized coordinate system. 



Arbitrary Constant Weak Fields 



For weak fields, we recognized in the previous sections that the assumption 
v ~ 1 or k 2 ~ 0 is consistent with the results for the velocity shifts in the low- 
frequency limit. To the same order of accuracy, we are now interested in the 
form of the polarization tensor for an arbitrary constant field configuration in 
the weak-field limit. As usual, the contribution proportional to IIq Pq V can 
be neglected. The physically more interesting part of which exists in 

the || JL subspace spanned by , can be written as a matrix: 



n ^ = 



n \\ 

0 \l v l v l ) 



(2.121) 



Here we have employed the fact that = z\/(a 2 + b 2 ) 2 + 0(k 2 ) > 0, and 
■y/ujj vj_ is the normalization of the tensor . Diagonalization of this matrix 
is straightforward; we obtain the eigenvalues 



A 1,2 




+ n± ^(i7||-i7x) 2 +402 V j|^ 



(2.122) 



With the aid of the weak-field approximations (2.110) to the functions IV, 
and setting = m 2 to our order of approximation, the functions II ||_j_ and 
0 can easily be evaluated: 



An = -2 



n ± = -2 



or Zk 



14 



m 4 a 2 + b 2 V 45 + 45 ^ 



or z k 



14, 



<9 = 2- 



m 4 a 2 + b 2 V 45 ' 45 

a 2 , 6 



ab ■ 



mr 



45' 



(2.123) 



Inserting these equations into (2.122), we find, surprisingly, that the eigen- 
values do not depend on the invariants a and b any more: 
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, a 2 /II , 3 

Ai,2 = -2 — j Zk — ± — 

m 4 V 45 45 



(2.124) 



Hence, we obtain exactly the same eigenvalues for the ||,_L subspace of the 
polarization tensor as in the crossed-fields case in the weak-field limit: II\ 2 = 
\i y 2 - Obviously, this leads to the same velocity shifts as obtained in (2.117): 



Vl,2 = 1 - 



11 _3_\ z ± 

45 45 / to 4 k 2 



(2.125) 



combinations of v and v^_: 



The eigenvectors u) 1 2 of the matrix (2.121) are constructed from certain linear 

(2.126) 



U 

u? o = m 



1,2 



* 1,2 



The coefficients can immediately be calculated: 



mi = 



V a 2 + b 2 
b 

V a 2 + b 2 



rrii = 



= 



V « 2 + b 2 ’ 

a 



Va 2 + b 2 

Inserting (2.127) into (2.126), the eigenvectors u^ 2 reduce to 



*Fk^ 



Fk» 



w7 = 



Un = 



yj %k yj Zk 

so the directions of the polarization modes are given by 

ei ~ — k X E + B (B ■ k)B , 
e 2 ~ k X B + E — (E ■ k)E. 



(2.127) 



(2.128) 



(2.129) 



This is the same result as the one that we found for the crossed-fields case. 
In the weak-field limit, the polarization tensor near the light cone does not 
depend on the complete set of invariants fc 2 , a, 6, Zk , but loses its dependence 
on the pure field invariants a and b. Furthermore, the remaining influence of 
the field via 

z k = k 2 [B 2 sin 2 6 B + E 2 sin 2 0 E + 2k • (B X E)] (2.130) 

is strongly dependent on the orientation of k, E and B. 



2.3 Induced Electromagnetic Current 

The calculation of the electromagnetic current induced by an external constant 
magnetic held is comprehensively performed in all technical details. For this, we 
solve explicitly the proper-time equations of motion, mapping Lorentz matrix com- 
putations onto simple Pauli spin algebra. The outcome agrees with Adler’s result, 
and coincides with the findings obtained via the polarization tensor in the previous 
section. 
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As an example of a proper-time calculation, we investigate here the electro- 
magnetic current induced by an external c-number field. For this, we employ 
the Hamiltonian formalism for a quantum mechanical one-particle problem 
as outlined in Sect. 2.1. In particular, we review the approach of Adler [4], 
whose calculation of the effects of photon splitting and vacuum birefringence 
is based on this technique. Adler’s work, in turn, was inspired by the ideas 
and calculations of Minguzzi [124]. A generalization of the results for vacuum 
birefringence was later worked out by Daugherty and Lerche [51]. 

As our physical starting point, consider the effective field equation of a 
wave field a M propagating in an external (constant magnetic) field: 

n2 <V = — O’m)- 

The appearance of the expectation value of the current on the right-hand side 
accounts for the fact that an external field modulated by a propagating wave 
field generally induces vacuum polarization, which can be diagrammatically 
understood in terms of virtual loop processes. Hence, the expectation value 
of the current depends, in general, nonlinearly on the external field as well 
as on the wave field. In a classical sense, the linear part with respect to the 
wave field can then be reinterpreted as an effective contribution to the index 
of refraction for the propagation of light. 

To be precise, if there is a contribution to the current which is linear in 

Un) - -/( w , B ) a to /( w > B ) arbitrary, 

then the field equation is solved by a plane wave field 

iw( — t + nfi-x) 

W /i C 1 



where the refractive index is given by 



= 1 - 



,.,2 



Bearing this concept in mind, we begin with the expression for the current 
given in (2.4): 



(0| j^(x) |0) A = ietr[y M G(x, x|A)]. 



(2.131) 



Representing the Green’s function as an operator in configuration space, 
G(x,x\A) = (x\("fll + m)~ l \x ), and treating the inverse operator expression 
symmetrically leads us to 



(0|^(:r)|0 } A =--tr 






m 2 — ( 7 U ) 2 m 2 — (7 n ) 2 



('yn)lx) , 

(2.132) 



where we have omitted two mass terms which contain odd numbers of 7 
matrices. With the aid of the proper-time integral representation (2.7) of the 
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inverse squared Dirac operator, we may write the current in the following 
way: 

OO 

(0|.7 /i (x) |0) A = C - Jds e~ lm2s tr [YY(A n »U (s) |x) +7 ! V i (x|77’(s)77„ M] > 
o 

(2.133) 

where the proper-time evolution operator U(s) was defined in (2.12). 
Introducing the proper-time definitions 

l*)= M 0 )): <®(s)| = (z(0)| U(s), 

x = x(0), 77 =77(0), (2.134) 

x(s) = U-\s ) x(0) U(s), 77(s) = U-\s) 77(0) U(s), 

and employing the identity j^Y = —g^’ 1 ' — ia llv (cf. Appendix A), we obtain 
our final representation of the current for arbitrary external electromagnetic 
fields: 

OO 

(J M (*))= - e -Jdse~ im2s tr [<x( S )|77^( S ) + 77^(0)|x(0)} 

0 

+icr MI/ (x(s)\II v (s) — II v (0) |x(0)}] . (2.135) 

To evaluate this expression further, one has to diagonalize the operators 
77(s) ± 77(0) in configuration space to obtain c-number expressions, i.e. one 
has to represent 77(s) ± 77(0) as functions of x(s) and x(0). This can be 
achieved by solving the one-particle equations of proper-time motion as given 
in (2.15). 

Obviously, analytical solutions to these equations for arbitrary electro- 
magnetic fields are difficult to obtain. Therefore, we specialize to the situ- 
ation under consideration in which a plane wave field / ^ propagates in a 
strong magnetic background field B pointing along the z axis. In particular, 
we make the replacements 



TV * -T //ii/(£), £ — TL • X, 



(2.136) 



where = (1, n) defines the direction of propagation of the plane wave and 



F^ = F^ = 



/0 0 0 0 \ 
0 0 50 
0 -B 0 0 
\0 0 00 / 



= i B <T t 2- 



(2.137) 



The index “t” (transverse) on the second Pauli spin matrix specifies the 
Lorentz subspace in which = 1,2. The remaining tensor space will 
be labeled by the index “1” (longitudinal). Introducing a potential a M = 
e At exp(io;^) for the plane wave field, with a polarization amplitude e^, we 
may write 
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f [iv — di/CL^ — iuj( 



- n u e 



V^[L 



)e 



io ;£ 



(2.138) 



An appropriate expression for 77(s) — 17(0) is immediately obtainable by 
integrating the equations of motion (2.15) once. Considering the second line 
of (2.15), we obtain 



d!7 At (s) 

ds 



= e 



{F[iv + f/j.v[£(s)}}II v 

+B l/ (s){F liV + / MV [^(s)]} + <9 M / A "[£(s)] 



= 2e n u (s)F /ll/ + ^(s), 



(2.139) 



where 

= e (n v i s ) //4£0)] + //4£0)] nv i s ) + A s ) d nf Xv [t,( s )^j ■ 

(2.140) 



Substituting the equation of motion (2.15) for £ M (s) into (2.139), we can 
perform a proper-time integration to find 







(2139) 





F, 



fils' 



Mt) 



= Cj-^[x„(s)-x v ( 0)] + y«P M (t)dt, (2.141) 

0 



where we have introduced the short form 1 v = eF^ . 

Some more effort is needed to find an expression for 77(s) + 77 (0) . Clearly, 
we have to extract the information from an integration of the proper-time 
equation of motion for the coordinate variable x For this, we need an inte- 
grate expression for the momentum 77 ^(s), which can be obtained by solving 
the momentum equation of motion (this time without inserting the coordi- 
nate equation of motion) : 



^^ = 2eC(-)^( S ) + ^( S ) 

7 7„(s) = (e 20 '"^ "77,(0) + J ( e 2 C<->( s -t)^ v $ v {t) dt. (2.142) 

o 



In the following, the 4x4 matrix calculations can be simplified by employing 
Pauli spin algebra, since it is useful to note that 



2C ( - } s (2 'i 37) 2i eBsa t2 

„2C ( -> 



= li + It cos 2eBs + icrt 2 sin 2eBs, 



(2.143) 
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where we have made use of the notation li = diag(l, 0, 0, 1) and l t = 
diag(0, 1, 1, 0). Owing to the absence of an electric component in the external 
field, the Lorentz tensor structure completely decomposes into purely trans- 
verse and longitudinal tensors which commute, i.e. [-PpQt] = 0 for arbitrary 
tensors P and Q. 

Substituting (2.143) into (2.142), we find the desired expression for TJ^s), 
which can now be inserted into the equation of motion for ^(s). Upon inte- 
gration, we obtain 



l (s)-x fl (0) = J dt = 2 J n^(t)dt 

o o 

s 

= 277^(0) J d t (li + l t cos 2 eBt + i(j t 2 sin 2eBi) /Ji 1 ' 
o 

s 

+2 J dr(li + l t cos2e5r + icr t 2 sin 2 eBr)^ x 
o 

r 

x J dt(li + lt cos 2eUt — icr t 2 sin 2eBt) x ''<P y (t). 



Integrating the r integral in the second term by parts leads us to 
eB[x^{s) - x M (0)] 

= [2eBs lj + l t sin 2 eBs + icr t 2 (1 — cos 2 eBs)] J II „ (0) 

S 

+ J dt [2eBsl\ + lt sin2eB(s—t) — ia t 2 cos 2 eB(s-t)] <P„(i) 
o 

s 

- J dt (2 eBs li - i 0 * 2 )/ $„(t). (2.144) 

o 



From this equation, we can determine 17^ ( 0) completely in terms of and 
tp ^ by multiplying by the inverse of the matrix 

M = 2 eBs li + It sin 2 eBs + iert 2 (1 — cos 2 eBs) . (2. 145) 



As can be easily checked by direct matrix multiplication, the inverse M 1 is 
given by 



...i 1 1 + cos 2 eBs „ i 

M = lj T It — — u”t2- 

2 eBs 2 sin 2 eBs 2 

Multiplication of (2.144) from the left yields 



(2.146) 



n x ( o) 

= (M^ 1 ) A /i eB[x M (s) - x M ( 0 )] 



(2.147) 
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+ dt 



- c it 



-ll - -It - icr t 2 
s 2 



1 + cos 2 eBs\ 
2 sin 2 eBs ) 



v Mt) 



ii + it 



(1 + cos 2 eBs) sin 2 eB(s-t) cos 2 eB(s-t) 



-10t2 



2 sin 2 eBs 2 

(l+cos2eJ3s) cos2 eB(s-t) sin2 eB(s-t) 
2 sin 2 eBs 2 






JA 



For the moment, let us go back to the integrated momentum equation of 
motion (2.142), use the representation (2.143) and add 77^(0) on both sides: 

n^ s ) + n,(o) 

= [21i + l t (1 + cos 2eBs) + icr t 2 sin 2eBs\ "77,, (0) (2.148) 

S 

+ j dt [li + l t cos2e7?(s — t) + iot 2 sin2eT?(s — t)]^"<£„(t). 
o 

Finally, we can insert the representation of 77„(0) found in (2.147) into the 
right-hand side of (2.148) to obtain 

S 

n„{s)+n^) = C^>[x v {s)-x v {dj\+ J dtT(s,t)^<P v (t), (2.149) 

0 

where, after a straightforward calculation, the tensors and T(s,f) are 
found to be given by 

C (+) = [21i + l t (1 + cos 2 eBs) + icr t 2 sin 2el3s] M” 1 eB , 

fv 0 0 0 ^ 

T = v li + C(s, t) l t + S(s, t) icr t 2 = 



. (2.150) 



0 C(s,t) S(s, t) 0 
0 —S(s, t) C{s, t ) 0 
\0 0 0 vj 

The functions in the last line of (2.150) are computed as 

2 1 _ sin eBvs . cos eBvs — cos eBs 

v= 1, C(s, t) = — — , S(a,t) = : .(2.151) 

s sm eBs sin eBs 

Now we are in a position to rewrite the current for this special field configura- 
tion. We insert the expressions for 77(s) ± 77(0) found in (2.141) and (2.149) 
into the representation for a general current, (2.135). Noticing that terms pro- 
portional to x(s) — x(0) vanish, since (x(s)|x„(s) — x„(0)|x(0)) = x v — x v = 0, 
we find the induced current to be 



O'/iO)) = 77 / ds e“ 



dt tr{ [T(s, t)/ + icr^^] (x(s) \$ v (t) |x(0)) } . 



(2.152) 
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It should be pointed out that this result is exact for an external field config- 
uration of a constant magnetic background field modulated by a plane wave 
field where both are treated without taking radiative corrections into account. 
Diagrammatically, we are dealing with an electron loop with infinitely many 
external photons of both origins (constant and plane wave fields) attached to 
it. 



However, c-number results can only be gained by performing further ap- 
proximations. The physical circumstances of a weak plane wave field propa- 
gating in a strong background field suggest expanding the current in powers 
of the plane wave amplitude. This can be achieved in a systematic way by 
going over to an interaction picture in which the “free” evolution operator 
contains only the constant magnetic field, and the plane wave field is treated 
as a perturbation. We therefore need the exact solution of the zeroth-order 
approximation, i.e. the case of a magnetic field alone - a problem that was 
originally solved by Schwinger [148]. 

In the following review of this solution, the zeroth-order quantities will 
be marked by a superscript In particular, we employ the definitions 

= H (0)2 - | aF, 

= U(0)-\s) 77 (o) (0)C/ (o) (s), 

= C/(0)- 1 (s)x (o) (0)C/ (o) (s), 

(2.153) 

= i-Ber t2 . 



UW(s) = e~ iHWs , 


#(°) 


77 (0) (0) = IJW, 


n(°\a) 


/ o 

III 

'o 


x(°\s) 


[x ( °\nl 0) ]= i g^, 




[n^,n^] = ieFpp, 


(F)^ 


It is evident that the equations of motion i 



daffi (s) 
d s 

d n { °\ s ) 

d s 



2 

2 eF^n^(s), 



(2.154) 



since [F Mi ,,77p°^] = 0. Solving these equations is a simple task of integration. 
In particular, the momentum equation yields 



nj?Xs)= (e 2eFs )/U( o) (0) 

= (li + l t cos 2 eBs + i cr t 2 sin 2 eBs)^ v II ^ (0) 
= i?(s)/ 17^(0), 

/I 0 0 0\ 

0 cos 2 eBs sin 2 eBs 0 
0 — sin 2 eBs cos 2 eBs 0 
\0 0 0 1 / 



where we have defined 

F(s)/=R(s) =e 2C(-) * = 



(2.155) 



(2.156) 
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Inserting (2.155) into the coordinate equation of motion (2.154), the latter 
can immediately be integrated, with the result 



c^(s)=x^(0) + 2(l lS + l t 



sin 2 eBs 



cos 2 eBs — 1 



2 eB 



icr t2 - 



2 eB 



= f(o)+/( s );j7(°)(o), 



o) 

(2.157) 



where 






= \{s) 

(2 s 
0 



(2.158) 



0 0 0 \ 

(sin 2 eBs)/(eB) (1 — cos2 eBs)/(eB) 0 
0 (cos2eBs — l)/(eB) (sin 2 eBs)/(eB) 0 

\ 0 0 0 2s J 

One further necessary ingredient for the following expansion is the zeroth- 
order evolution operator in coordinate space (x'°) (s )|x^(0)), which has al- 
ready been found in (2.39) of Sect. 2.1. In the case of a purely magnetic field, 
the expression reduces to 

(x( 0) (s)|x (0) (0)) = (x (0) | U (0) (s) |x (0) ) 

= e i(e/2)lTFs . (2.159) 

(47ts) 2 sin eBs 

This completes our investigation of the zeroth order, i.e. the constant-field 
part of the problem. 

Our next intermediate step is characterized by finding the relations be- 
tween the zeroth-order and exact quantities. While the exact and zeroth-order 
coordinates at zero proper time coincide and thus can be viewed as a bound- 
ary condition, the corresponding canonical momenta differ by the plane wave 
amplitude: 

*m(0) = 4 0) (°)> ^m(O) = Pm - e \ = n l 0) (°) - erf 0) (0)], (2-160) 

whereby = n ■ x^. By insertion of (2.160), the exact Hamiltonian can 
be related to the zeroth-order Hamiltonian: 



H= H 2 - - a(F + f) = {iZ(°>( 0) - ea M [^°)(0)]} --a(F+ f) 






a^°\o)}n^^ + n^a^°\o)\ 



(2.161) 



- e a / ,[e (o) (0)]a^K(°)(0)] 



' [IV 



n (0) m 



In order to divide the complete evolution operator into the “free” (zeroth- 
order) and interaction parts, 

,-i H^s ; 



U(s) = t/ (0) (s)t/i(s) 
we employ the operator identity 



e~ iHs = e 



'CAM, 



(2.162) 
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~A+B — A 



At T3„At 



= e A T exp / d te At Be 



(2.163) 



where A = — i H^s and B = — i (77 — H^^s. Obviously, the interaction 
evolution operator can be computed as 

Ui(s) = T exp ^-i J d tU(° r \st) (H - H^)sU {0 \st) 

= T exp ^-i J duU^~\u){H - ijW)C/ (0) (u)J . (2.164) 

With the help of (2.161), the integrand in the exponent of (2.164) can be 
written as 



-i £/ (0) \u)(H - H^)U^(u) 

o#.[f (0) («)] ^ (0)/ » + 77^(0) a^ 0 >(u)] 



= ie 



(2.165) 



-ea M [£<°>(u)K[£<°>(u)] + ^ ^)/H£ (0) K>] 



where we have repeatedly inserted factors of 1 (m). Incidentally, 

the proper-time-dependent er matrix in the last term is given by 

0$(u) = uM~\u) U^\u) = Q~ i i e / 2 '> ucrF e ^/^ F , (2.166) 

since [77^ 0) , cr MJ/ ] = 0. 

Substituting (2.165) into (2.164) produces 



U\ = T exp< ie j du 

^ o 



a M [£< 0 >(u)] n^(u) + 77^(0) a M [e (0) («)] (2-167) 



^ea^ 0 >(u)K(^(u)) + ^{u)r\^\u)} 



The collection of the tools required for the interaction picture is now com- 
pleted, and we can turn to the remaining problem of rephrasing the ob- 
ject (x(s)|^ y (t)|x(0)), which is contained in the equation for the the current 
(2.152): 

(x(s)|<7v(7)|x(0)) 

= (x(0)\U(s)U~ 1 (t)^>^(0)U(t)\x(0)) 

= (a: (0) | £/(°) (s) f/i (s) Uf 1 (t) {7^°)“ 1 {t)<&v (0) (t) Ui (t) | a:(0) ) 

= (ar<°> (s) | t/i (s) C/r 1 (t)^S 0) (*) (*) | a:(0) <°> > . (2.168) 
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We may give a more explicit expression for (t) in the last line by employ- 
ing its definition in (2.140) and noting that II 1 ' = — ea v \ 

^\t) = e(nM*{t) f, p [^°\t)} + f„ p [e°\t)} nW»{t) 

—2 ea'[tW(t)]U P [Z W (t)] + | CT <°>(t)^/ A -K (0) (t)]) .(2.169) 

Inserting the last line back into (2.168), which finally enters (2.152), we obtain 
the desired expression for the current in the interaction picture: 



Unix)) 



OO S / 

- c — J dse~ lm2s J dftr < [T(s, t)/ + itr /i I/ ] 

n n ^ 



(2.170) 



x (a:^ (s) | f/i (s) Uf 1 (t) 



n i 0 )p (t) U P [^°\t)] + uie o) (t)] n Wp (t) 



-2 ea p [^ 0 \t))f vp [e°\t)] + 



C7i(t)|* (o) (0)) 



To study the birefringence properties of the vacuum, it is sufficient to consider 
the first-order approximation of the current, which describes the lowest-order 
interaction of the plane wave field with the magnetically induced currents. 
Since the operator in the large square brackets in (2.170) is at least of first 
order in the plane wave field, the interaction evolution operators U\ have to 
be approximated simply by U\ ~ 1. Furthermore, the third term in the large 
square brackets of (2.170) has to be omitted, because it is already of second 
order in the plane wave field. The resulting first-order approximation of the 
current then reads 



OO s , 

(Ux)) = -y J ds e~' m2s J dt tr j [T(a, f)/ + i<] (2.171) 



x(x(s)| 



n p {t)U P m\ 



+um\ n p {t) + 



i*m. 



For convenience, we have omitted the superscript in (2.171) and shall 
continue to suppress it in the following, since it should be clear that from 
now on we are only dealing with zeroth-order quantities. 

Our final task is to evaluate the matrix elements of the operators in the 
large square brackets in (2.171) in coordinate representation and then take 
the Dirac trace to obtain an explicit c-number expression. With the aid of 
(2.155), and inverting (2.157), we obtain a coordinate representation for the 
momentum operator: 
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n(t) = R (t) ■ r\s) • [x(s) - a?(0)], 



(2.172) 



where we have used matrix notation and the dot denotes matrix multiplica- 
tion. Similarly, we can represent the plane wave phase by 



(2.157i 

exp[iw£(t)] = exp[iwn • x(t)\ '= exp{iw n ■ [x(0) + 1(f) • 77(0)]} 



= exp|A[x(s)] +S[x(0)]|, 
where we have defined 

A[x(s)] = i ton - 1(f) • l” 1 (s) • x(s), 

S[x(0)] = icon ■ [l — 1(f) ■ I - 1 (s)] • x(0). 

Since the commutator of x(s) and x(0) is nonvanishing, 

[x(s) M , ac(0)^] (2 '= 57) -i/ F (s), 
we find for the commutator of A and B 
R(s,t) = -^[A[x(s)),B[x(0)]} 



(2.173) 

(2.174) 

(2.175) 



-^[n-\(t)] ■ |n- [1 - 1(f) ■ I \s)]} 



(2.176) 



which is c-number valued. Employing the Baker-Campbell-Hausdorff for- 
mula, 

e A+B = e A e B e~^ A,B \ 
we may write for (2.173) 

gi «>£(t) _ e 4[i(s)] e B[i(o)] e “ 2 fl(s,t) i (2.177) 

Regarding the products of 77 p (t) and in (2.171), we need to order 

the operators in such a way that all the x(s)’s stand to the left of all the 
x(0)’s. In this way, they can be pulled out of the (x(s)| and |x(0)) brackets. 
For this, we make use of the operator identity 

e a b = be a + [a, b] e a , 

which we apply to 



x“(0) e A [ x ( s )] = g^K®)] x“(0) -co\n- l(f)]VI lW l 



(2.178) 



and 



e B > l( 0 ) l^)“ = i a ( S )e fl > l(0)1 

+cj ^l(s) • |n • [1(f) • l _1 (s) — l] e B [ x(0) l. (2.179) 

With the help of (2.172)-(2.179), we can rewrite the equation for the current 
(2.171): 
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OO S s 

0'/*( a: ))=-y w Jdse~ lm2s Jdte^^^tr 1 (T/ + icr/) (x(s)|x(0)) 

0 0 ^ 

Vr(I) • jr^s) • [n- !(*)]+«• \(t) ■ r^s) -n}^ (2.180) 



+ 2 ri c/ AK<7 (^« 



where the transition amplitude (x(s)|x(0)) = K(s) is given in (2.159). Here 
it should be pointed out that the object cr(t)\ K can only be pulled out of 
the sandwiching bracket to the right, since (x(s)| possesses a nontrivial Dirac 
structure, while |x(0)) does not. 

In principle, the Dirac traces, as well as the matrix manipulations, could 
be performed in full generality, but for reasons of lucidity we shall make use 
of the simplifying specifications of our particular configuration. The magnetic 
background field points along the positive 2 axis; we consider the plane wave 
as propagating along the positive x direction: = (1, 1, 0, 0). Consequently, 

we denote the two possible polarization vectors as e(J) = (0,0, 1,0) and ey 1 = 
(0, 0, 0, 1), specifying the directions perpendicular (_L) and parallel (||) to the 
plane spanned by the magnetic field and the propagation direction. 6 This 
choice also normalizes the plane wave potential, = a| ±IJj ± . 

The matrix multiplications can be performed straightforwardly; l(s) is 
given in (2.158), and its inverse reads 



1 , eB sin2eBs 

— — lj T — -It 

2s 2(1 -cos 2 eBs) 



-eH ot2, 



(2.181) 



which enables us to evaluate the function R(s,t) defined in (2.176): 



R(s, t) 



^(1-u 2 ) 



cos eBvs — cos eBs 
eB sin eBs 



(2.182) 



Using (2.156), the plane- wave-independent part of the first term in the second 
line of (2.180) is now computed as 



R(f)-jl 1 (s) • \n ■ 1(f)] +n- 1(f) 





/ ~ v \ 


\ Q — 


C(s,t) 


r~ 


S(s,t) 




\ 0 / 



~ ( n ' T ) c 



(2.183) 



where the functions S(s,t) and C(s,t ) are defined in (2.151) and the matrix 
T is given in (2.150). 



Besides the use of different metric conventions, note that Adler’s definition of 
perpendicular and parallel components is opposite to ours, since he refers to the 
magnetic field vector of the photon field. In the present work, the polarization 
vector designates the direction of the electric field vector of the photon field. 
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Substituting (2.183) and (2.159) into (2.180) leads us to 



'd s 



27 r / s 2 



a 

= i 



e -im 2 S / dteW 2 fl( S ,t) 



1 eBs 
4 sin eBs 



(2.184) 



x t [T ■ f-(n ■ T)] tr U e/2)s<TF | +i[f ■ (n-T)] tr 



e 



Vfj(e/2)s<rF 



+ ^(T-n)/ AK tr 



:h i(e/2)saF 



<?\ K (t) 



~2 tr 



CTu e 



! 'e i (e/2)so'-F 



O' A K (l) 



The simplest way to, on the one hand, calculate the Dirac traces and, on 
the other hand, extract the desired information is to project the current onto 
the polarization directions by multiplying by aj) Incidentally, the current 
can only be constructed as a Lorentz vector from the polarization vectors, 
because gauge invariance implies transversality. 

For the reduced case of a purely magnetic field, the exponential function 
of crF can be conveniently represented by Pauli spin matrices (cf. Appendix 
C): 

e i(e/2) s <rF = cos e Bs + i cr 3 sin eBs, (2.185) 



where it is understood that the (2 x 2) Pauli matrices are doubled blockwise 
in Dirac space. We can immediately read off that 



tr 



g i(e/2 )saF 



4 cos eBs , 



(2.186) 



and, in conjunction with the identities listed in Appendix C, we also find 



tr 



n 0 i(e/ 2 ) SCTF 



—4 sin eBs 



(2.187) 



The same relation holds for the third trace in (2.184). For the last trace, we 
observe that, by employing the contraction of the current with aj*j ± ^, this 
object can be written as 



tr 



(flU a ) 



J(e/2)(s-t)a p 



(fl± a ) 



= — w 2 (±16),, ± (cos eBvs 



e i(e/2)tcrF 

- cos eBs), 



(2.188) 



which involves the definition of a^(t) in (2.166). (Details can be looked up 
in Appendix C.) 

After inserting these traces into (2.184) and computing the scalar products 
by projecting onto the polarization directions, we can summarize our findings 
for the induced current regarding this particular field configuration as 




-u} 2 A l± (u,B) a^ ± , 



(2.189) 



2.3 Induced Electromagnetic Current 



45 



where Aim is given by 



= SLj ^e im2s J d te“ 2R ^ % X (M), (2.190) 

o o 



and the functions JVim are constructed from simple combinations of trigono- 
metric functions: 



„ T eBs cos eBvs „ , 9 sm e.Bus\ 

Nu = eBs cot eBs I 1 — v z + v ) 

" sm eBs \ sin eBs ) 



N ± = - 



eBs cos eBvs veBs sin eBvs cot eBs 



sin eBs sin eBs 

2e£?s(cos eBvs — cos eBs ) 
sin 3 eBs 



(2.191) 



(2.192) 



An expression for R(s,t) has been obtained in (2.182). 

These results are completely consistent with the findings of the previous 
section, in which the same situation was discussed via the polarization tensor. 
To be precise, the induced current as calculated here for a plane wave field 
moving in an external constant magnetic field is related to the polarization 
tensor given in (2.91) according to 



0’|pj_) = ~n> iV {ujn) a v (wn), (2.193) 

which simply reflects the implications of the fundamental principle of lin- 
ear response as described by (2.2). At this point, we want to stress that 
the present calculation has been simplified in two ways compared with the 
previous section: first, we assumed a specific propagation of the plane wave 
field, namely, perpendicular to the B field; however, the general case can be 
immediately restored by inserting a factor of sin 2 6 into the function R(s,t ) 
in (2.182) and as an overall factor in front of (2.190) ( 0 denotes the angle 
between the magnetic field and the propagation direction). Secondly, we as- 
sumed during the whole calculation that the plane wave field propagates on 
the light cone by requiring that n M is a null vector (cf. (2.136)). This implies 
that the refractive index for the plane wave field remains close to 1, which 
is, of course, not true in the strong-field limit (B/B CI > n /a). As has been 
demonstrated by Daugherty and Lerche [51], it is not essential to maintain 
this assumption during the calculation; indeed, from the details given above, 
it is obvious that it only enters the final matrix manipulations, which sim- 
plify slightly for this case. Of course, the strong-field limit has to be analyzed 
without any a priori assumptions about the wave vector (cf. Sect. 3.2). 

As mentioned above, the technique outlined here can also be applied to 
the calculation of photon-splitting matrix elements [ ] . For this, the induced 
current has to be evaluated to second order in the plane wave field. 

Finally, let us mention that a generalization of the present calculation to 
the case of parallel electric and magnetic fields was treated by Daugherty and 
Lerche [51]. 
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2.4 Schwinger’s Equivalence Theorem 
and the Axial- Vector Anomaly 

Schwinger’s equivalence theorem, which claims equivalence between the pseu- 
doscalar and pseudovector interactions of a pseudoscalar field with electrically 
charged fermions at low energies, is derived. It is shown that there are no non- 
perturbative effects arising from the all-order coupling between the external elec- 
tromagnetic field and the fermion loop. 

The relation between the equivalence theorem and the axial-vector anomaly is 
outlined, proving their coincidence for small photon momenta only. As a corollary, 
we find that there are no nonperturbative contributions to the axial- vector anomaly 
to one loop in the domain of small photon momenta. 

Besides being a powerful tool for investigating QED, proper-time techniques 
provide for a detailed study of axial couplings between pseudoscalar and 
fermionic fields in the presence of external electromagnetic fields. In particu- 
lar, problems related to gauge invariance can be investigated thoroughly. 

In his seminal work, Schwinger [148] proved the equivalence theorem , 
which states that, in the low-energy regime, a pseudoscalar interaction be- 
tween a spinless neutral meson and a fermion field leads to the same result 
for the decay of the meson into two photons as a pseudovector interaction. 

In the following, we shall reinvestigate Schwinger’s equivalence theorem 
with special emphasis on the nonperturbative domain, i.e. the coupling of the 
external electromagnetic field to the fermion loop to all orders. Furthermore, 
we shall point out in what respect the equivalence theorem is related to the 
axial- vector anomaly as discovered by Adler, Bell and Jackiw [2, 2 ]. 

2.4.1 Equivalence Theorem 

Let us first describe Schwinger’s view on the two-photon decay of the neutral 
pion, which was inspired by a paper by Steinberger [156] on the question of 
the equivalence of various interaction Lagrangians. In fact, it will turn out 
that the present considerations cannot be immediately applied to pion decay, 
owing to a mismatch of energy-momentum regimes; let us, however, stick to 
the historical viewpoint for the time being and later translate our findings 
into modern language. 

Pseudoscalar Interaction 

For the pseudoscalar interaction between a neutral meson field and a fermionic 
field, we begin with the Lagrangian: 



where g denotes a dimensionless coupling constant, and (j> represents the 
spinless meson (pion). The fermion was identified with the proton; nowadays, 
ip should be associated with a quark appearing in three colors. 




(2.194) 
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In order to describe the decay of the pion into two photons, Schwinger 
replaced the fermion fields by their vacuum expectation value in the presence 
of an external electromagnetic field: 

£ PS £el = -i gHx) \{[^{x),l^{x)]) A 

= g 4>(x) tr j 5 G(x, x\A), (2.195) 

where taking the vacuum expectation value in the first line translates into a 
time-ordering prescription via the point splitting procedure: 

itr7 5 G(a;, x\A) = lim 7 5 (0|T^(a; , )r)>(a;)|0) A . (2.196) 

x'-^x 

In (2.196), the limit has to be performed with respect to the time- like com- 
ponents, and the invariant distance ( x — x') 2 should be space-like, i.e. > 0. 
Equation (2.195) can be diagrammatically represented as in Fig. 2.2, 




Fig. 2.2. Diagrammatic representation of (2.195) 



where the double line represents the coupling to the external electromagnetic 
field to all orders. Obviously, taking the vacuum expectation value with re- 
spect to the external field corresponds to integrating out the fermions to one 
loop. 

Inserting the proper-time representation of the Greens function G (cf. 
(2.7)) into (2.195), we obtain 



( OO 

75 (x\(M — 7 17 )i j d S e- is [ M2 -^ n ) 2 ]| 
o 

( °r . 

= gM (j){x) itr I 75 / dse -lsM (x|e ls ^ |x) 



(2.197) 



where M denotes the mass of the fermion. To arrive at (2.197), we have 
employed the fact that traces of odd numbers of 7^’s together with a 75 
vanish. In the last term of (2.197), we encounter the trace of the proper-time 
evolution operator. 
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Now Schwinger, identifying the internal fermions with protons, argued 
that the momentum of the outgoing photons of the pion decay is much smaller 
than the mass of the loop particle. Therefore, the electromagnetic fields as- 
sociated with the photons vary slowly compared with the length scale set 
by the Compton wavelength of the proton. As a consequence, the constant- 
field approximation for the proper-time transition amplitude appears to be 
appropriate in the present situation. 

However, it is not the proton whose fluctuations constitute the loop; 
rather, the loop particles have to be identified with the quarks. 7 But the 
current masses of the light quarks are much smaller than the momenta of the 
outgoing photons of the decaying pion, which completely spoils the slowly- 
varying-field assumption. The fact that pion decay is nevertheless describable 
with our final results can be attributed to the general form of the anomaly, 
which is already revealed by the constant-field approximation. This will be 
elucidated in more detail later on. 

Hence, let us forget about pion decay and simply proceed with the 
constant-field/low-photon-energy approximation assuming that a heavy fer- 
mion runs in the loop. With the aid of (2.39), we insert the representation of 
the evolution operator into (2.197) and find 8 






(4tt) 2 



o 



eas ebs 

~ sin eas sinh ebs 



x tr{ 7 5 e i(e/2),TFs } 

oo 

= -H-0ji<Kx)e 2 J <ise-““ 2 ab 

= --±4,( x )g = -±4,(x)E-B 
7 T M 7 T M 



1« _9_ 
4 7T M 



ct>{x)F^F^. 



(2.198) 



This equation represents the low-energy effective Lagrangian of a pseu- 
doscalar interaction between a spinless meson and a heavy fermion in an 
external field. Although we have included the coupling of the loop fermion 
to the external field to all orders, the final result is only of second order 
in the electromagnetic field strength. Hence, if we had expanded the loop 
perturbatively in a, then only the graph with two external photons would 
have contributed to the final result. Note also that we have encountered no 

7 In fact, in order to obtain a reasonable value for the pion lifetime, one has to 
take the number of colors into account. 

8 Remember Q = (l/^F^F^ = ~E • B = ab. 
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singular terms while calculating G{x,x\A)\ the dangerous terms vanished by 
Dirac algebraic arguments. 

Pseudovector Interaction 

Let us consider a pseudovector interaction of the following particular form: 

£PV = _i 2 iff Yi ■ (2.199) 

Classically, this pseudovector interaction Lagrangian is formally equivalent 
to the pseudoscalar counterpart as defined in (2.194), since 

£ PV ' = P ' + [V>,75#]} +s.t., 

e '= m ' -i 9<t>(x) i[^(x),7 5 t/)(x)] + s.t., 

where “s.t.” means “surface terms”. Here, we have first integrated by parts 
(i.b.p.); then we have employed the equations of motion (e.o.m.). However, 
at the quantum level, things become more complicated. Proceeding in the 
same way as in the pseudoscalar case, we naively arrive at 

£eff = d^ix) ^ ([^(x), 757^(20] ) A 

= -i^|y^</>(x)tr757' i G(x,a:|H) 

1 = P ' i-^^(x)d li tT'ys^G(x,x\A) +s.t. (2.200) 

Now we are in trouble! Not only do we have to face the problem of singularities 
in G(x,x\A), but we also have to give a meaning to the derivative at this 
singular coincidence point. The first guess is to introduce an appropriate 
point splitting regarding the two arguments of the Green’s function. But 
then, we have to keep in mind that G{x,x'\A) is a gauge variant quantity 
while G(x, x\ A) is not. In order to ensure gauge invariance, we have to replace 
the ordinary derivative with covariant derivatives. Following Schwinger, we 
reinterpret the critical term in the last line of (2.200) as 

d/xtr 757 ^ G(x,x\A) (2.201) 

-► x }^ x {\^ ~ ieA n( x ')] + K + ie A^{x")\ }tr 757^ G{x' ,x"\A). 

One can easily check that the right-hand side reduces to the left-hand side 
after naively taking the limit in a formal sense. Now we could follow Schwinger 
and evaluate the right-hand side of (2.201) in the weak-field limit, i.e. up to 
second order in the field strength. This would again correspond to a triangle 
graph, which is known to contribute solely to the axial-vector anomaly to 
any finite order of perturbation theory. 
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Instead, we shall try to maintain the coupling to the external field to all 
orders in order to pursue the question of possible nonperturbative contribu- 
tions to the meson-photon interaction. Of course, the price we have to pay 
is that we are strictly tied to the slowly-varying-field approximation. 

So, let us employ the representation of the fermionic Green’s function in 
an arbitrary constant electromagnetic field given in (2.46): 



G(x, x'\A) 

OO 

= ${x,x'\ A )j^ j ^ rn- + et]^(x-x’Y 



x exp 



— i M 2 s — L(s ) + -(x— x')f(s)(x— x') 



exp ^i-erFs^ , (2.202) 



where 



f(s) = eF coth(eFs), 






eas ebs 



sin eas sinh ebs 



and 



<P(x, x'| A) = exp 



X 

ie J d 4 (a»® + - x%) 



(2.203) 



completely carries the gauge dependence of the propagator. Having separated 
the gauge dependence in this way, we may also write 



G(x' , x" | A) = <2>(x' , x" | A) G(x' , x" | A sf ) , 



(2.204) 



where G(x', x"|Asf) is the Green’s function evaluated in the Schwinger-Fock 
gauge and depends only on the field strength. 

Regarding the required Dirac traces, we find 

tr 757 m e i ( e/2 )< jFs = 0 (odd number of y^'s) , 

tr 757^7 « e i(e/2)<TFs = 4i(eT*F> ia + P g» a - eT* j , (2.205) 

where we have extensively employed the results of Appendix C. Furthermore, 
it is useful to find an explicit expression for the function f(s) a p; for this, we 
make use of the spectral representation of F as described in Appendix B. 
We obtain 

f(s) a p = a 2 1 b 2 ( fl2 9<*P + F ap) ebcothebs 

H — 2 , , n ( b 2 g a 0 - Flp) ea cot eas. (2.206) 

a z + b z H 

Note that f(s) a p is symmetric. In view of (2.201), we also need the result of 
the following derivative construction: 
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{ K - ieA M (x')] + [d” + i eA^(x")\ }$(x',x"\A) 
(2.203) . e + Ii^ A (y _ ieA A1 (x / ) 

-ie ^4 M (x") + i j d£ K F% j + ieA M (x") 
= ieF^x' - x") A ^(a; , ,x"|A). 



<5(0/,®" | A) 



(2.207) 



It is this term, in particular, which we would not have discovered in the naive 
calculation of the left-hand side of (2.201). 

Let us now combine our findings from (2.201)-(2.207), leading first to 

{ [% - i eA^{x')\ + [<9" + i eA^x")] }g(x',x"\A) 

(2=04)(|[^_ieA M (x / )] + [S"+ ieH^x")]}^', x"|A)) G{x' , x"|A SF ) 

+^( x ', x "\ A )^ + c>")g(x',x"|Hsf). (2.208) 



Since G(x' ,x"\Asf) = G(x' — x"|Asf) owing to translational invariance, the 
second term vanishes completely in the constant-field case (we shall come 
back to this point later on). Upon insertion of (2.205) and (2.207) into (2.208), 
we arrive at 



tr 757 ^{ — ieA^(x')] + [5" + ieA^x")] }g(x', x"\A) 



le 



2(4tt) : 



;F^{x'-x"Y${x',x"\A) 



(2.209) 



-je lM2s (/ + eF) a p(x' — x"Y 



i ebs 



x exp 



r (x'-x")f(a/-x") 



sin eas sinh ebs 
4i (eT *F t * ot + P g^ a — eT* F ^ . 



Here, we encounter the Lorentz product, 

~(x' - x") A F X/ ,feT*F^ a + P g^ a — eT* F^ (/ + eF) a0 (x' - x"f. 



Since this product is symmetrically contracted with {x' — x"), only symmetric 
terms with respect to the field strength indices can contribute. With the aid 
of the explicit representation of the function / (2.206), we find, after a tedious 
but straightforward calculation, 



-{x' - x") x F Xll (eT*F^ a + Pg» a - eT* i^*“)(/ + eF) a0 (x' 

( sin eas sinhe6s N 



/ / //\A 

(x — x ) 



b 2 



i 2 b 2 



\ sinh ebs sin eas 



9 A/3 



x'Y 

( 2 . 210 ) 
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. sin eas 2 s i n h e bs 
sinh ebs sin eas 



F- 



A/3 



0 x'-x'Y ■ 



Now, the crucial observation is that the result of (2.210), read together with 
the last factor of the second line of (2.209), can be related to the total deriva- 
tive of the function f(s) a p with respect to the proper time: 



d f(s) a0 (2.206) 
ds 



a 2 + b 2 



2 t2 

a b 



sinh 2 ebs sin 2 eas 



9x/3 



1 



1 



K 3 



. , , , • , , ,, ■ ( 2 . 211 ) 

sinh ebs sin easy 

Inserting (2.210) and (2.211) into (2.209), we obtain a comparably simple 
representation of the desired expression on the right-hand side of (2.201): 

tr 757 M {[^-i eA„(x')\ + [<9" + i eA^{x")\^G{x' ,x"\A) 



— ^(x 7 , x"\A) ab dse lM s exp 



1 { f //\r/ t "\ 

-(x -x )f(x -x ) 



x(x'-x")“ - 



(- - df ^ s ap ) {x'-x"f (2.212) 



= i—*F^F^$(x',x"\A) I dse~ iM2s ^-exp 
2n / ds 



i( x' — x")f(x' — x") 



Substituting this result back into the starting point, i.e. into the effective 
Lagrangian in (2.200), yields 

C p J=-l^^xyF^ (2.213) 



x lim < 

x' ,x" — >x 


f oo 

<2>(x',x"| A) fdse~ iM2s ^-exp 


\{x'-x") f(x'-x") 


1 




, 0 




J 



Comparing this with our result for the pseudoscalar interaction in (2.198), it 
is obvious that an equivalence exists between the two different interactions 
on the quantum level if the limit expression in (2.213) finally reduces to 1 for 
any kind of constant electromagnetic field. By construction, the proper-time 
integration has to be performed before the limit x ' , x" — > x can be taken. For 
example, if we interchanged this processes in (2.213), then we would find a 
zero result, since (d/ds)e° = 0. 

Employing the fact that <P(x r , x"\A) — > 1 for x',x" — > x, we are left with 
the question of whether 



lim [ dse~ iM2s ^-exp 
x',x"^x J ds 

0 




x")f(x' 
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Lb = P ' lim U- iM ' s 

x' ,x" — >x 



exp 



-( x '- x 'y( x '- x ") 



J ) 0 



+i M 2 lim / dse” iM2s exp 



X' .x" — >x 



-(x'-x")f(x'-x") 



(2.214) 



equals 1. First, we shall demonstrate that the boundary term vanishes. For 
this, note that convergence at the upper bound s — > oo is ensured by the 
prescription M 2 — > M 2 — ie. This prescription can also be implemented by 
rotating the contour of s slightly below the real axis. Hence, the boundary 
term at s — > oo vanishes because the mass term approaches zero exponen- 
tially. For the lower bound, we need an expansion of f{s) a p for small values 
of s; referring to (2.206), we obtain 

f(s) a p = - g a @ + — s F 2 p + 0(s 3 ). (2.215) 

In the limit s — > 0, only the first term is relevant, and we find that 
exp [(i/4)(x' — x") 2 / s\ — > 0 for s slightly below the real axis and (x' — x") 2 > 
0 space-like. This completes the proof that the boundary terms in (2.214) 
vanish. 

Hence, we are left with the following integral: 9 



I{a 1 b) = \M 2 lim 



OC 

J 



dse- iM2s exp 



4 <*'- 



x")f(x' 




(2.216) 



To get a better feeling for this integral, let us first make contact with 
Schwinger’s original work and consider the weak- field limit. Since the result- 
ing effective Lagrangian is already of second order in the fields (cf. (2.213)), 
it suffices to consider the integral (2.216) in the zero-field limit; for this we 
employ the expansion (2.215), because the weak-field expansion of f{s) a p 
coincides with the small-s expansion: 



OO 

I (a = 0,6 = 0) = M 2 lim / dse _iM2s exp 

x' ,x" — >x J 



i (x' -x") 21 



s — ► — IS 



= M lim dse~ ms exp 
x' ,x" — >x J 



4 s 
1 (x' - x") 2 



^ lim M (x' — x") Ki [M (x' — x")l 

x' ,x" — >x 
= 1, 



(2.217) 



9 It is amusing to see that a (forbidden) interchange of the limit and the s inte- 
gration would lead to the correct result: (2.216) — * 1. But this is accidental, as 
we shall soon demonstrate. 
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where we have encountered the modified Bessel function K\(x ), which has a 
simple pole at x — > 0 with residue 1. Inserting this into (2.213), we obtain 
the effective Lagrangian of a pseudovector interaction between a spinless 
meson and a heavy-fermion field in the presence of a slowly varying and 
weak external electromagnetic field: 

£ ° PffV = Hx) F ^ v - (2 - 218) 

This is identical to the outcome for the pseudoscalar interaction and shows 
the essence of Schwinger’s equivalence theorem for the low-energy regime. In 
this sense, the terminology “low energy” refers to the energy of the outgoing 
photons (variation of the field strength) as well as the strength of the field. 

But we want to go one step further and prove the validity of the equiva- 
lence theorem without the weak-field assumption. For this, we have to show 
that I(a, b) = 1 for all values of a and b. As stated above, interchanging the 
limit and the s integration in (2.216) is not a valid operation. To put this 
in mathematical language, let ft( x) be the limiting value of a sequence of 
functions ft n (x); then 



O O 

/ h(x) dx = lim / h n (x) dx 

n *°° J 



if the h n {x) are integrable on the interval M = [a, b ], and 
lim || ft - h n || M = 0, 

n — ►oo 

where 



(2.219) 



\\h\\ M = sup {| ft (a:) | | x € M} 



is the supremum norm of h(x) for the complete interval M = [ a,b ]. As an 
example, let us study the zero-field limit, as treated above; here, ft(s) = e -lM s 
and ft(a;'_a:")(s) =e -lM s ed/ i )i x ~ x ) / s . Obviously, we find for the supremum 
norm: 



lim 



„-i M 2 s 


„-i M 2 s „ 


[i Or'-x") 2 ] 


e 


— e exp 


4 s 



M 



= 1 , 



( 2 . 220 ) 



where M denotes the s integration interval from zero to infinity slightly 
below the real axis. This does not satisfy the necessary criterion (2.219) for 
an interchange of limit and s integration. But since a direct evaluation of the 
integral in I(a,b) in (2.216) is not in sight, we have to find an indirect way 
to prove that I (a, b) = 1. 

For this, we choose to work in the special Lorentz frame where the electric 
and magnetic field are antiparallel, say B = ae z , E = —be z . In this case, 
the square of the field strength tensor appearing in (2.206) has the form 
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F 2 af} = 



( -b 2 0 0 0 \ 

0 —a 2 0 0 

0 0 —a 2 0 ' 

\ 0 0 0 b 2 ) 

In this special system, the exponent ( x ' — x")f (x' — x") takes the particularly 
simple form 

(V — x") a f(s) a p(x' — x"Y = Ax| e&coth ebs + Ax(j_ ea cot eas, (2.221) 

where 



Axjj 1 = ({x' — x")° , 0, 0, (x r — a:") 3 ) ; 

Ax± = (0, (x' — x") 1 , (x' — x") 2 ,0). (2.222) 

As an exercise, let us briefly investigate the case of purely magnetic fields, 
i.e. b — > 0. Then, we may write 



(x 1 -x”) a f(s) a/ 3 (x' -x")^ 



= Axn — [-Ax, ea cot eas 
6=0 11 s 



(2.223) 



= Ax 2 



- Ax i 



ea cot eas — 



where we have introduced the short form Ax = x' — x" . For /(a, &), we obtain 
in this limit 



l(a,0) s ~= ls lim M 2 dse 
Ax 2 Ax 2 ^o J 



-M 2 s 



(2.224) 



x exp 



1 Ax 2 
4 s 



— - Ax^ ( ea coth eas — 



Here, we have replaced the limiting process x',x" — > x, which is essentially 
four limiting processes, one for each component, by two limiting processes 
Ax(j_, Ax 2 — ► 0. Now, we make the crucial observation that the Ax^ — > 0 limit 
can be interchanged with the integral, since 



lim 

Az^o 



1 — exp 




^ea coth eas 




(2.225) 



Equation (2.225) is a consequence of the fact that eacoth eas— 1/s is bounded. 
Therefore, we obtain, following (2.217), 



/(a,0) 



OO 



lim M 2 
Ax 2 ^0 



! 



dse 



— M 



exp 




= 1. 



(2.226) 



Hence, for arbitrarily strong external magnetic fields, the equivalence theorem 
remains valid. But (2.226) tells us more; since a function which is constant 
(= 1) on the whole positive real axis is constant on the whole complex plane, 
we obtain 
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1 = /(«.()) n _. ; = I(—ib, 0) = 7(0,6), (2.227) 

which can be easily checked with the aid of the definition of I(a, b) in (2.216). 
Thus we do not even find a further contribution to the low-energy effective 
Lagrangian of the pseudovector theory for arbitrarily strong electric fields. 
The latter might have been expected, since the present situation resembles the 
familiar situation of calculating the effective Heisenberg-Euler Lagrangian for 
constant electric fields, which is known to reveal a nonperturbative imaginary 
part. 10 

The final proof that /(a, b) equals 1 is a simple generalization of the pre- 
ceding considerations. For this, we write for (2.216), with the aid of (2.221), 



/(a, 6) 



OO 

lim iM 2 [ dse~ iM2s exp ( 
Ax 2 ,Ax^.,Ax?,— ►() J \ 



i Ax 2 
4 s 



(2.228) 



x exp 



i 2 / 

— Accii e&cothe&s — 



exp 



Ax, ea cot eas — 



Let us now treat the Ax 2 1 Ax 2 L and Ax 2 limits independently of each other, 
i.e. let us forget about the fact that Ax 2 = Ax\ + Ax 2 . Note once again that 
the convergence of (2.228) is ensured by the prescription M 2 — > M 2 — ie. 
Furthermore, assuming that Ax 2 7 Ax 2 ±: Ax 2 , > 0, we are allowed to rotate the 
s integral contour by an angle of — 7r/4 below the real axis: 



— e i7r / 4 s = e -i7r / 4 s = — t =(1 — i) s. 

V2 



The convergence properties are maintained, since 
e“ 

i Ax 2 ' 



- iM2s ' exp ( --^=(1 + \)M 2 s ) , 



exp 



4 s 



exp - 



1 (1 — i) Ax 2 



(2.229) 



(2.230) 



Now the crucial observation is that both combinations, eb coth ebs — 1/s and 
eacot eas — 1/s, are bounded along the new s path. Hence, we are allowed to 
interchange the Axy and Axx limits with the s integration, leading to 



/(a, b) = lim ie' 71 ’/ 4 M 2 / dse~ ie W4m2s exp 

Ate 2 — >0 J 

o 



i7r/4 



Ax 2 



= lim M 2 [ dse~ M2s exp(~ — ) 

Ax 2 ^o J \ 4 s ) 

o 

10 It can in fact be checked directly that 1(0,6) has a vanishing imaginary part by 
summing all the residues of the cot poles on the lower imaginary s axis. One 
finds that the residues are all of the order of Ax| and thus vanish for Ax| — > 0. 
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= 1. (2.231) 

This proves Schwinger’s equivalence theorem for arbitrary constant external 
field strengths. Besides, it proves that there are no nonperturbative contribu- 
tions to the effective pseudovector Lagrangian in an external field. This latter 
remark will be important for the discussion of the axial-vector anomaly be- 
low; it demonstrates that there are no nonperturbative corrections to the 
anomaly at the one-loop level in the low-energy domain. To summarize the 
equivalence theorem, the main result can be written in the form 

£ e p i = 4? it m F ^ F " V = £PffV ’ (2 - 232) 

where the tree-level interaction Lagrangians for the pseudoscalar and pseu- 
dovector interactions are defined in (2.194) and (2.199). Equation(2.232) 
holds for arbitrary electromagnetic field strengths as long as the fields vary 
slowly compared with the Compton wavelength of the fermionic loop particle. 

2.4.2 Axial- Vector Anomaly 

It has often been emphasized in the original literature [2, 84, 107] that the 
discovery of the axial- vector anomaly (Adler-Bell-Jackiw anomaly) has its 
roots in Schwinger’s work, which we have outlined above. Hence, we want to 
investigate the connection between the equivalence theorem and the axial- 
vector anomaly more closely. The anomaly can be summarized in the state- 
ment that the axial vector current is not conserved, not only because of an 
explicit breaking of axial symmetry by mass terms, but also because of the 
appearance of the F llv *F tlv term induced by quantum effects. 

In order to derive this statement, we have to give a proper definition of 
the axial- vector current, as well as its divergence. On the classical (or even 
on the operator) level, the following definitions are reasonable: 



35 = ^1/0,75^], 


(2.233) 


35 = \ [0. 757 m 0] ■ 


(2.234) 




(2.235) 



The investigation of the pseudoscalar interaction leads us to the conclusion 
that for the expectation value of the axial (scalar) current, we obtain 

(j 5 ) A = itr 75 G(*,*|A) = - i ^±F^F^, 

^ £e! = -i 9 (js) A ■ (2.236) 

Note that {j$) A is well defined, since the singularities of G(x,x\A) are re- 
moved from {j§) A by the Dirac algebra. 

Applying the same ideas to the case of the pseudovector interaction, we 
are tempted to write 
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“( 3s) A = itrqsq ' 7 G(x,x|A)". (2.237) 

However, this time, the singularities of G(x,x|A) are not protected by the 
Dirac trace, which causes (j£) A of (2.237) to be an ill-defined quantity. 

Since we are not interested in (j£) A but in its divergence, we shall not pur- 
sue the question of how (2.237) can be corrected; instead, let us concentrate 
on a proper definition of d^{j^) A . 

During the study of the pseudovector interaction, we were forced to give 
a proper meaning to tr 757 ^ G(x, x\A), which then was achieved by a 
point-splitting procedure in (2.201). Separating the Green’s function into 
its gauge-dependent part <P(x' . x"\A) and a part which depends only on the 
field strength G(x',x"|Asf) (cf. (2.204)), we arrived at (2.208): 

x , [d'^ - ieA^x')} + [<9" + ieA„{x")\ }g(x', x"\A) 

(2.208) lim j-j e p x (x' _ x") A $(x', x"\A)\ G(x', x"| A S f) 

x' ,x" — >x 

+ lim $(a!y\A)(ff„ + &^G{rf,a!'\A SE ). (2.238) 

We remark that the second term is singularity free, since the singularity in 
G(x' , x"\Aq,y) must be of the form ]im. x > !X n- t . x G(x 7 , x"|Asf) ~ lim x -' )X "_, x 
h(x' — x"), with some function h, where h{ 0) — > 00 . Hence, the singularity 
vanishes in the second term because of the symmetry in d ' <-> d”. Since 
( P(x' , x"\A) — > 1 for x',x" — > x, we are led to interpret the second term in 
(2.238) as a building block for an appropriate definition of the divergence of 
the axial- vector current: 

dn(j^(x)) = lim ^(x 7 , x"\A) (9^+9") itr 757 ^ G(x 7 , x ,, |H S f)-(2.239) 

We want to stress that this definition is gauge-invariant, although it is not 
obvious, since, after taking the limit x' ,x" — > x, the left-hand side and the 
first term on the right-hand side of (2.238) are separately gauge-invariant. As 
shown above, the trace of 75 7 ^ times the first term on the right-hand side of 
(2.238) yields 

lim ieF^x (x 7 - x") A ^(x / , x"\ A) tr 757 ^ G(x 7 , x"\ A S f) 

x' ,x" — >x 

= ^ . (2.240) 

Z7T 

Substituting (2.239) and (2.240) into the right-hand side of (2.238) and apply- 
ing the equivalence theorem to the left-hand side, i.e. tr [qsq^LHS of (2.238)] 
= —2M(j§) A , we finally arrive at 

dvUg) = -2i M(j 5 ) + F^r (2.241) 

This is the well-known equation for the divergence of the axial- vector current 
and represents the Adler-Bell-Jackiw anomaly. However, by employing the 
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equivalence theorem, we have proved (2.241) only for the special case of 
constant external fields. Hence, Schwinger’s work on the constant-field case 
is only capable of deriving the anomaly in a certain energetic regime, namely 
the low-energy domain. 

Within the usual diagrammatic approach [2], one also arrives at (2.241), 
but this time without assuming that the electromagnetic field has to be con- 
stant. For example, starting with the left-hand side of (2.240), one may ex- 
pand the Green’s function in an arbitrary external field as a power series in 
the coupling constant: 

G= S F + ieS F ~fAG 

= Sp + ie Sf 'yA Sf + (ie) 2 Sf 'yA Sf 'yA Sf + . . . . (2.242) 

Then one observes that the zeroth-order term vanishes owing to the Dirac 
algebra, the second term leads to the correct anomaly and the e 2 term, as 
well as all higher-order terms, vanishes in the limit x',x" — > x, since these 
terms do not contribute to the singular behavior of G. (Note that the latter 
statement holds only perturbatively.) It is exactly this second term in (2.242) 
which is diagrammatically represented by the famous triangle graph. 

To summarize our considerations, we want to stress that the equivalence 
theorem corresponds to the axial-vector anomaly in the domain of slowly 
varying fields only. This is the only regime where the equivalence theorem is 
applicable. Moreover, by giving up the weak-field assumption, we were able 
to prove the equivalence theorem for arbitrary field strength, i.e. to all orders 
in the external field. Translating this into the language of the anomaly, we 
have proved that there are no nonperturbative corrections to the anomaly in 
the case of slowly varying external fields. 

Without the constant-field approximation, one can show perturbatively 
that (2.241) for the divergence of the axial- vector current holds for arbitrary 
fields to any finite order of perturbation theory. The price that one has to 
pay for this stronger statement is that one loses control over the nonper- 
turbative domain. Within these techniques, the question still remains open 
as to whether the anomaly receives nonperturbative corrections for rapidly 
varying fields. However, Fujikawa’s path integral procedure, which is truly 
nonperturbative for arbitrary A fields [81], demonstrates that (2.241) is not 
modified. 

Let us remark in this context that the validity of (2.241) for arbitrary 
electromagnetic fields does not tell us anything about the question of whether 
the equivalence theorem holds for arbitrary photon energies. 

Let us conclude this chapter with an interesting observation for the 
constant-field case. Inserting our findings for (js)" 4 for constant fields (cf. 
(2.236)) into (2.241), we find that the divergence of the axial- vector current 
vanishes: 

d,U0 = o. 



(2.243) 
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This is already clear from its definition (2.239), since G(x' , £"|Asf) for con- 
stant fields depends only on x 1 — x" . Hence (2.239) vanishes simply by an- 
tisymmetry. Nevertheless, the result appears rather unfamiliar, because it 
signals the conservation of the axial- vector current at the quantum level, al- 
though this current is not conserved at the classical level owing to the explicit 
breaking of the axial symmetry by mass terms. Therefore, the constant-field 
case is an exceptional situation which creates an “inverse anomaly”: a clas- 
sically and explicitly broken symmetry is restored by quantum effects. 

Since the fields are considered to vary slowly compared with the Compton 
wavelength of the fermionic loop particle, one can reinterpret the constant- 
field limit as a heavy-fermion limit. In this sense, very heavy fermions, al- 
though breaking the axial symmetry very violently, do not violate the axial- 
vector current conservation after quantum fluctuations are taken into ac- 
count. 

Finally, let us stress once more that our considerations remain strictly at 
the one-loop level. Similarly to the Fujikawa method, we do not take photonic 
fluctuations into account; hence our statements about nonperturbative contri- 
butions to the anomaly do not touch questions about higher-loop corrections. 
In particular, we cannot comment on the correctness of the Adler-Bardeen 
theorem [3]. 



3. Nonlinear Electrodynamics: 
Effective- Action Approach 



In the preceding chapter, we studied examples of quantum electrodynamic 
effects induced by external fields by means of the electron Green’s function, 
the polarization tensor and the induced current. The “operator” language we 
employed made explicit reference to the quantum nature of the underlying 
physics. In contrast, the physical effects that we extracted from this formal- 
ism, such as optical birefringence, resemble closely the features of classical 
physics. 

An important tool that interpolates between the full quantum theory and 
classical field theory is provided by the concept of the effective action: once the 
high-energy degrees of freedom which are invisible in the classical low-energy 
domain are integrated out, the resulting effective action for the remaining 
degrees of freedom can be employed to define a new quasi-classical theory. 1 
In particular, classical equations of motion can be derived whose solution 
provides for the description of quantum effects in the language of classical 
physics. In this sense, electrodynamics represents an excellent laboratory for 
investigating and applying this concept. 

While classical electrodynamics is distinguished by Lorentz invariance, 
gauge invariance and linearity, quantum electrodynamics features only the 
first two principles. Nonlinear effects clearly are of quantum origin and arise 
as a result of the polarization of the vacuum. Self-interactions of the electro- 
magnetic field are a by-product of the creation and subsequent annihilation 
of charged pairs. Hence, these self-interactions are nonlocal in nature on a 
microscopic scale, but are, of course, causal, since the interaction with the 
charged pairs is causal. Nevertheless, this nonlocal nature can be disregarded 
on a scale where the virtual creation and annihilation processes can no longer 
be resolved. This scale is set by the Compton wavelength of the electron, 
A = 1/m. 

Taking these considerations into account, the one-loop effective Lagrang- 
ian of QED for constant fields as derived in Sect. 2.1, (2.49), read together 
with the classical Maxwell Lagrangian, can be employed to define a local 
classical field theory that enriches classical electrodynamics by nonlinear self- 

1 In this work, we do not elucidate the question of whether the effective action 
might be used as a defining action for a new quantum theory of the residual 
degrees of freedom. This type of approach is reviewed in [ 11] and [ 36]. A 
subtle example demonstrating the limits of this approach can be found in [32]. 
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interactions of the fields. The validity of this theory is, on the one hand, 
limited by the validity of the loop expansion; on the other hand, it is restricted 
to scales larger than the Compton wavelength of the electron; in particular, 
the fields must vary only slowly on that scale. 



3.1 A First Look at Light Propagation 



It is proposed that the study of light propagation in a perturbed QED vacuum is a 
useful tool for investigating the manifold properties of such a quantum ground state. 
As an introductory example, we solve the equations of motion for low-frequency 
light propagating in an external field for the simplest nonlinear extension of classical 
electrodynamics: the Heisenberg-Euler Lagrangian. We obtain a modified light cone 
condition, and discuss polarization dependences as well as velocity shifts. 

The main purpose of this section is to illustrate on a simple level the various 
techniques and necessary approximations that will become important in the sections 
to follow. 

Since we are interested in the borderline at which classical and quantum 
physics meet, i.e. where nonlinear effects begin to play a significant role, the 
study of the propagation of light, which is trivial in the classical vacuum the- 
ory, serves as an important example. In this introductory section, we briefly 
present our strategy by means of the simplest example: the lowest-order cor- 
rection to the classical Lagrangian. The more sophisticated investigations 
in the later sections are based on the same philosophy as described in the 
following. 

In (2.50), we found the lowest-order correction to the Maxwell Lagrangian; 
in combination, the “next-to-leading-order” classical electrodynamics can be 
defined by the lowest-order Heisenberg-Euler Lagrangian 



C = —J- + ci T 2 + C 2 G 2 , 



8 a 2 

Cl 45 m 4 ’ 



_ 14 or 
° 2 45 m 4 



The equations of motion derived a la Euler-Lagrange reduce to 
dC „ dC „ dC dF a0 _ dC 



0 = 



dA t 



d>1 d(dM d, *dF af3 d(dM 



= -2d L 



dF u 



(3.1) 



(3.2) 



where the differentiation with respect to F must be performed under the 
constraint of antisymmetry of the indices. With the aid of the identities 



^ = i pa p _9G_ = K pa0 

dF a p 2 ’ dF a0 2 

we end up with 

0 = <9 m ( F ^ - 2d T F^ - 2 c 2 Q *F tlv ) , 



(3.3) 

(3.4) 



which represents the vacuum field equation for the lowest-order nonlinear 
electrodynamics. This equation is accompanied by the usual Bianchi identity 
df/F 1 - 1 '" = 0 . 
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Applying these equations to the case of propagating light in a background 
field, we decompose the field strength tensor into the background part F ^ 
and the plane wave field f^ v , F — > F ^ . We assume the background 

field to be constant or, at most, slowly varying and disregard any derivative of 
F^ v . Furthermore, we linearize the field equations with respect to the plane 
wave field; this is equivalent to neglecting self-interactions of the propagating 
light itself. In a certain sense, the plane wave is thus reminiscent of a test 
charge in classical electrodynamics. It is important to note that the lineariza- 
tion does not necessarily imply that the background field is much stronger 
than the plane wave field. The zero-field limit for the background field will 
be well defined. 

With these assumptions, we obtain for the field equations, with the aid 
of the Bianchi identity, 

o = d^r - d F a 0 F^ d^f a0 - c 2 *F a0 *F ^ d^r 0 . (3.5) 



Since f^ v is assumed to be a plane wave field, we perform a Fourier 
transformation to momentum space, where can be written as = 
k^a v — k^a^. Here, = (w,k) denotes the wave vector of the propagating 
light; consequently, the wave vector replaces the derivative in (3.5); d M — > ifc'L 
Finally, we impose the Lorentz gauge constraint on the gauge potential a M of 
the plane wave, = 0, and arrive at the light cone condition 

0 = k 2 a" - 2 ci F a 0 F^ u k^k a a 0 - 2 c 2 *F a0 *F ^ k^k a a 0 . (3.6) 

Obviously, we recover the trivial light cone condition k 2 = 0 in the classical 
limit where Cy 2 = 0. For the two different polarization states that solve this 
field equation, we may try the ansatz 

< ~ W = *F^ x k x , ~ Fk» = F^kx, (3.7) 



which is motivated by our findings in Sect. 2.2, (2.128). Employing the fun- 
damental algebraic identity (B.4), we find upon insertion of this ansatz into 
(3.6) 

: 0 = (k 2 - 2 c 2 z k + 2 c 2 Tk 2 )*Fk v - (2 Cl Gk 2 )Fk v , 

«£ : 0 = (k 2 - 2ci z k )Fk u - (2c 2 Gk 2 )*Fk v , (3.8) 

where z k = (F^k ^) 2 denotes the Lorentz- and gauge-invariant quantity that 
we have already encountered in (2.63). Owing to the last terms on the right- 
hand sides of (3.8), it is obvious that the polarization states a ± and a% do 
not solve the field equation exactly, however, since we are in the weak-field 
domain, we expect k 2 to deviate only slightly from the trivial light cone 
condition: k 2 = 0 + 0(ci i2 ). Hence, the terms proportional to Ci^Qk 2 or 

c-iTk 2 will be of higher order in Cp 2 and should be omitted for reasons of 

consistency. We find 



cii : k 2 = 2c 2 z k , 



k 2 = 2ci z k - 



(3.9) 
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Introducing the phase velocity v via k ^ = (/c°, k) 
we obtain 



a 



n . 
1 • 



a 



m . 
2 • 



vi 



V2 



i _ £fc _ 1 _ if 

2 k 2 45 m 4 k 2 ’ 

1 _ i £fc 

1 k 2 45 m 4 k 2 ’ 



k(w,k), where k = jk|, 



(3.10) 



which exactly corresponds to our findings in (2.125). An appropriate repre- 
sentation of the invariant Zk in terms of E and B fields is given by (2.130): 
z k = k 2 [B 2 sin 2 0 B + E 2 sin 2 0 E + 2k ■ (B X E)\ . 

Within the Lagrangian formalism developed so far, this result was first 
obtained by Bialynicka-Birula and Bialynicki-Birula [28] (correcting and gen- 
eralizing a result of Klein and Nigam [115]). For purely magnetic fields, a sim- 
ilar calculation was performed by Adler [ ]. With a noncovariant notation, 
the weak-field computation is nicely reviewed in [125]. 

The following remarks should be made: 



(1) The light cone condition derived from the field equations is of second or- 
der in k ^ by construction, since the effective Lagrangian (3.1) contains no 
derivatives of the field strength. The latter property, in turn, follows from 
disregarding the nonlocalities of the effective action. As a consequence, 
the phase and group velocities of the propagating plane wave coincide, 
because the phase velocity does not depend on the frequency. This agrees 
with the fact that the results are only valid for soft photons, i.e. for low 
frequencies, oj <C m, which cannot resolve virtual processes. 

(2) The phase velocities do not exceed the vacuum velocity c (=1), because 
the coefficients c \ i2 are strictly positive. However, even if they were not, 
we could not draw the conclusion that causality is violated, because the 
signal velocity is equivalent to the front velocity; but the latter is related 
to the phase velocity in the infinite-frequency limit w — > oo (cf. the 
discussion of causality in Sect. 3.3.4). 

(3) From (3.8), we could read off the linear combinations of Fk M and 
which, as the correct polarization states, solve the field equations exactly 
[28]. In the following, however, we shall be particularly interested in light 
cone conditions averaged over polarization (and propagation direction), 
which might be called “sum rules”. In this case, the calculation following 
(3.6) can be simplified by multiplying (3.6) by a„ and summing over the 
two polarization states according to the rule 

(3.11) 

pol. 

This yields 

0 = 2 k 2 (l -(- 2c 2 E'j 2(ci + c 2 ).£fc. 

Neglecting the higher-order terms in ci j2 again, we finally find 



(3.12) 



3.2 Light Cone Condition 
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fc 2 = (c i+c 2 ) Zk =► t; = l-i(c 1 +c 2 )g-l-ii^ T g.(3.13) 
In fact, this is the arithmetic average over the polarization states of (3.10). 



3.2 Light Cone Condition 

We derive the light cone condition for low-frequency light propagating in an ar- 
bitrary constant external field. The dynamics of the field is determined by a La- 
grangian which is assumed to depend on the two invariants of the electromagnetic 
field but is otherwise left arbitrary. We show that the deformation of the light 
cone is determined by the (vacuum expectation value of the) energy-momentum 
tensor of the field, and a field-dependent factor which we call the effective-action 
charge. By inserting the one-loop QED effective action, the light cone condition 
is applied to the cases of weak fields and strong and superstrong magnetic fields, 
yielding refractive indices for the various electromagnetically perturbed vacua. We 
observe that a vacuum modified by a superstrong magnetic field shows striking 
similarities to a magnetized plasma. This domain is also studied with the aid of a 
renormalization-group-improved effective action. 

The comparably simple procedure of the previous section will be generalized 
in the present section in order to obtain a light cone condition for a plane 
wave propagating in an electromagnetic field of arbitrary strength. Again, the 
basic quantity will be the effective action of QED in an arbitrary constant 
background field. The success of this strategy will further support the picture 
of considering the modified vacuum as a medium. 

In the following, we mimic the steps of the preceding calculation, but this 
time drop the weak-field assumption for the background field. To be concrete, 
we first demand that the following condition is satisfied: 

(I) Any length scale which we can construct from the parameters of the system 
(except from the field strength) is large compared with the Compton wave- 
length of the lightest massive particle (the electron) A c = 1/m. This ensures 
that the effective Lagrangian is local, implies that the field is “slowly” varying 
and restricts the plane wave to be of “low” frequency (compared with X c ). 

This property is sufficient for considering the effective Lagrangian to be 
a function only of the gauge and Lorentz invariants of the electromagnetic 
field; a convenient choice for the two linearly independent invariants is given 
by the standard (pseudo) scalars 

T = l -F lw Fi JI ' = i(£ 2 - E 2 ), Q = l -F lw *F“' = -E B. (3.14) 

For example, the Maxwell Lagrangian can simply be written as Cm = — T . 
The general Lagrangian therefore reads 

C = C(F,G). (3.15) 

Note that a parity-conserving fundamental theory (e.g. standard QED) de- 
mands that C is an even function of Q. 
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Following (3.2) and (3.3), we obtain the equations of motion from C by 
variation: 

r)C BC 

0 = - ax; = a ^ c ^ + a ‘ c ( 3 . 16 ) 

where d p, dg denote the partial derivatives with respect to the field strength 
invariants (3.14) (and should not be confused with space-time derivatives 
9 m ). If we take advantage of the Bianchi identity and move d M to the right, 
we arrive at 

0 = (djrC) , (3.17) 

where MZ is given by 

= F^F a0 (d$rC) + *F a0 (d 2 gC) 

+d?gC (F^ *F a0 + *F^F a0 ) . (3.18) 

Note that M is antisymmetric in the upper, as well as the lower, indices: 

= -K$ = M Z~ 

Next, we apply this field equation to the case of a plane wave field 
propagating in a background field F^ v , i.e. we decompose the field strength 
according to F ^ — * F^ + In doing so, we specify the second basic 
condition as follows: 

(II) Vacuum modifications (= self -interactions) of the plane wave field are 
negligible; this is formally achieved by a linearization of the field equation 
with respect to f . 

Furthermore, neglecting derivatives of the background field (owing to con- 
dition I) , we can replace the partial derivatives by the wave vector of the plane 
wave field in momentum space; this leads us to 

0 = {d^C)k,r . (3.19) 

As mentioned above, we are primarily interested in an average over polariza- 
tion states. Representing the plane wave field strength by f^ = k^a 11 — Fa'' 
(imposing the Lorentz gauge = 0), multiplying (3.19) by a„ and sum- 
ming over polarization states according to (3.11), therefore yields 

0 = 2 (d F C) k 2 + MZ k^k 01 . (3.20) 

Equation (3.20) already represents the desired light cone condition and actu- 
ally indicates that the familiar k 2 = 0 will, in general, not hold for arbitrary 
Lagrangians. Our final task is to put MZ into a convenient shape. Using the 
powerful fundamental algebraic relations of the field strength tensors given 
in Appendix B, (B.3) and (B.4), we obtain 

MZ = 2 [(1/2 )F^F au (d% + d 2 g )C + 5*{Qd Fg C - Fd 2 g C)} . (3.21) 



